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MULTIPLY-REFINED ENUMERATION OF 
ALTERNATING SIGN MATRICES 

ROGER E. BEHREND 

Abstract. Four natural boundary statistics and two natural bulk statistics are considered 

for alternating sign matrices (ASMs). Specifically, these statistics are the positions of the I's 

f~^ ■ in the first and last row and column of an ASM, and the numbers of generalized inversions 

CN , and —I's in an ASM. Previously- known results for the exact enumeration of ASMs with 

Q_i' prescribed values of some of these statistics are reviewed in detail. A quadratic relation which 

^ ■ fully determines the generating function associated with all six statistics is then obtained. 

The derivation of the relation involves combining the Desnanot-Jacobi determinant identity 

with the Izergin-Korepin formula for the partition function of the six-vertex model with 

domain-wall boundary conditions. 
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1. Introduction 

A major focus of attention throughout the history of alternating sign matrices (ASMs) 
has simply been the derivation of results related to their exact enumeration. Such results 
typically state that the number of ASMs which satisfy specific conditions, such as having 
prescribed values of certain statistics or being invariant under certain symmetry operations, 
is given by an explicit formula or generating function, is equal to the number of combinatorial 
objects of some other variety satisfying specific conditions, or is equal to a number which 
arises from a particular physical model. 

A few examples of such results, with references to conjectures and initial proofs, are as 
follows: a formula for the total number of ASMs of arbitrary fixed size but with no further 
conditions applied (Mills, Robbins and Rumsey [102, 103], Zeilberger [145], and Kuper- 
berg [92]); a formula for the number of ASMs with a prescribed boundary row or column 
(Mills, Robbins and Rumsey [102, 103], and Zeilberger [1 IG]); formulae for numbers of ASMs 
invariant under certain natural symmetry operations (Robbins [123, 121], Kuperberg [!l^>], 
Okada [111], and Razumov and Stroganov [119, 120]); equalities between numbers of certain 
ASMs and numbers of certain totally symmetric self-complementary plane partitions (Mills, 
Robbins and Rumsey [104], and Fonseca and Zinn- Justin [To]); equalities between numbers of 
certain ASMs and numbers of certain descending plane partitions (Mills, Robbins and Rum- 
sey [103], and Behrend, Di Francesco and Zinn- Justin ['"]); equalities between numbers of 
certain ASMs and numbers associated with a certain case of the 0(1) loop model (Razumov 
and Stroganov [116], and Cantini and Sportiello [31]). For reviews of some of these results, 
see, for example, Bressoud [22, 24], Bressoud and Propp [25], Hone [78], Zeilberger [1 IS], or 
Zinn- Justin [149]. 

The main result in this paper is of the type which provides an expression for the generating 
function associated with numbers of ASMs of arbitrary fixed size with prescribed values of 
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particular statistics. In any ASM, the first and last row and column each contain a single 
entry of 1, with all of their other entries being O's, so that the positions of these I's provide 
four natural statistics which describe the boundary configuration of an ASM. The generating 
function under primary consideration in this paper is associated with these four boundary 
statistics, together with two statistics which depend on the bulk structure of an ASM, namely 
the number of generalized inversions (which will be defined in (3)) and the number of —1 
entries. Elementary properties of this generating function and these statistics will be outlined 
in Section 2. 

All previously-known results for the exact enumeration of ASMs have involved fewer than 
six of these statistics. A review of such results will be given in Section 3. That section 
will aim to provide a comprehensive account of all previous appearances of these statistics 
in exact enumeration, together with an outline of some further related results which can be 
obtained straightforwardly from the previously-known results. Recently, a result involving all 
four of the boundary statistics, but neither of the bulk statistics, was obtained independently 
by Ayyer and Romik [12, Thm. 2]. That result will be reviewed in Section 3.3. 

In ASM enumerations, it is common to refer to a certain order of refinement which is based 
on the number of boundary statistics involved. Hence, the primary generating functions used 
in this paper, and by Ayyer and Romik [12], can be described as quadruply- refined. 

The main result of this paper. Theorem 1, will be stated in Section 4, and proved in 
Section 5. The approach used in the proof is essentially as follows. First, certain standard 
techniques, which have played a crucial role in proofs of many of the other known enumerative 
results for ASMs, are used to obtain a determinantal expression related to ASM generating 
functions. More specifically, a bijection between ASMs and configurations of the statistical 
mechanical six-vertex, or square ice, model with domain-wall boundary conditions (DWBC) 
is used to derive a relation between ASM generating functions and the partition function 
of the model, and the Izergin-Korepin formula is used to provide an explicit determinantal 
expression for this partition function. Next, the Desnanot-Jacobi determinant identity is 
applied to the matrix in the Izergin-Korepin formula, which leads to a quadratic relation 
involving the required quadruply-refined ASM generating function, and generating functions 
associated with no ASM boundary statistics or with two ASM boundary statistics, for all four 
pairs of adjacent boundaries. It will also be seen that this relation, together with simple initial 
conditions, fully determines the six-statistic quadruply-refined ASM generating function. 

The key step in this derivation is the use of the Desnanot-Jacobi identity together with 
the Izergin-Korepin formula, and the reason that this leads to a relation of the previously- 
described form can be understood relatively easily. The Izergin-Korepin formula involves a 
certain prefactor, which depends on so-called spectral parameters mi, . . . , m„ and vi, . . . , f„, 
and the determinant of a certain n x n matrix, whose entry in row i and column j depends 
on Ui and Vj. Furthermore, in terms of generating functions for nxn ASMs, the parameters Ui 
and Vj are, in a certain sense, associated with row i and column j of the ASMs. For 
an n X n matrix M, the Desnanot-Jacobi identity (in the form used here) consists of a 
relation involving det M and five minors of M: the four connected [n — 1) x [n — 1) minors. 
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corresponding to deletion of the first or last row and first or last column of M, and the central 
(n — 2) X (n — 2) minor, corresponding to deletion of the first and last row and column of M. 
Therefore, it seems that applying the Desnanot-Jacobi identity to the matrix in the Izergin- 
Korepin formula, and using an assignment of spectral parameters in which -Ui, «„, Vi and f„ 
remain arbitrary, while U2, ■ ■ ■ , Un-i are equal and t;2, . . . , Vn-i are equal, should lead to an 
expression involving a quadruply-refined generating function for n x n ASMs, four doubly- 
refined generating functions for (n — 1) x (n — l) ASMs, and an unrefined generating function 
for (n — 2) X (n — 2) ASMs. Indeed, this essentially does occur, with some further, subsidiary 
aspects of the resulting expression being related to the form of the prefactor in the Izergin- 
Korepin formula, and to the fact that each corner entry of an ASM is associated with two 
of the parameters Ui, Un, Vi or t>„. 

A similar method, involving the combination of a different form of the Desnanot-Jacobi 
identity with the Izergin-Korepin formula, was used by Behrend, Di Francesco and Zinn- 
Justin [ ], as part of the proof of a result concerning a doubly-refined enumeration of ASMs 
and descending plane partitions. 

The Desnanot-Jacobi identity has, in fact, been associated with ASMs since these matri- 
ces first arose. In particular, the identity is used in Dodgson's condensation algorithm [60] 
for determinant evaluation, and it was through studies of this algorithm by Mills, Robbins 
and Rumsey that ASMs initially appeared. In particular, it was shown by Robbins and 
Rumsey [125, Eq. (27)] that if Dodgson's algorithm is applied to an n x n matrix M using 
a modified form of the Desnanot-Jacobi identity containing a parameter A, then the result- 
ing so-called A-determinant of M can be expressed naturally as a sum over n x n ASMs. 
The latter expression involves the two ASM bulk statistics which are used in this paper, 
and, as will be shown in Section 3.2, the A-determinant shares certain features with the 
quadruply-refined ASM generating function studied in this paper. For example, the modi- 
fied Desnanot-Jacobi identity satisfied by the A-determinant is analogous to the quadratic 
relation satisfied by the quadruply-refined ASM generating function. For further informa- 
tion regarding Dodgson condensation, A-determinants, and related matters, see, for example, 
Abeles [1], Bressoud [22, Sec. 3.5], Bressoud and Propp [25], Propp [11 1], or Robbins and 
Rumsey [125]. 

Finally, it should be noted that the Desnanot-Jacobi identity, the A-modified Desnanot- 
Jacobi identity, and further related equations studied by Robbins and Rumsey [125], and 
others, are closely connected with the Laurent phenomenon, cluster algebras, and associated 
areas. For reviews of such connections, see, for example. Hone [7b, Sec. 5], or Propp [ i... 
Sec. 10]. 



2. Definitions and Basic Properties 

2.1. Statistics. In this section, the standard definitions of ASMs and certain statistics for 
ASMs are given, and some simple properties of these matrices and statistics are identified. 
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An ASM, as first defined by Mills, Robbins and Rumsey [102, 103], is simply a square 
matrix in which each entry is 0, 1 or —1, and along each row and column the nonzero entries 
alternate in sign and have a sum of 1. 

It follows that, for any ASM A, 

Each partial row sum ^2]'=! ^ij' ^^^ each partial column sum Yll'=i ^i'j is or 1. (1) 

It can also be seen that any permutation matrix is an ASM, and that, in any ASM, the first 
and last row and column each contain a single 1, with all of their other entries being O's. 

For each positive integer n, denote the set of all n x n ASMs as ASM(n). For example, 
for n = 1,2, 3, these sets are 

ASM(1) = {(1)}, 
ASM(2) = 



(2) 




ASM(3) = 

For any A G ASM(?7,), define statistics which depend on the bulk structure of A as 
z/(A) = y^ Aij Aiiji, /^(A) = number of — I's in A, 

l<i<i'<n 
l<j'<j<n 



(3) 



and define statistics which describe the configuration of A at its top, right, bottom and left 
boundaries as, respectively, 

Pt{A) = number of O's to the left of the 1 in the top row of A, 

Pr{A) = number of O's below the 1 in the right-most column of A, 

Pb{A) = number of O's to the right of the 1 in the bottom row of A, 

PhiA) = number of O's above the 1 in the left-most column of A. (4) 

The statistics of (4) are given diagrammatically by 



t 

PL (A) 



<-Pt(^) ■* 

/O 1 0\ 





1 
A 1 




VO 1 0/ 



t 

PR (A) 



(5) 
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The statistic z/(A) in (3) is a nonnegative integer for any A G ASM(n), since it can be 
written as i^{A) = X]ri=i(Si'^i ^j'i)(Sj'=i ^ijO) ^^^ (1) implies that each factor in the 
summand is or 1. This statistic can also be written as i^{A) = J2i<i<i'<n- i<j'<j<n^ij ^i'j'i 
where this can be obtained from the definition in (3) using the fact that each complete row 
and column sum of A is constant. As will be seen in Section 5.2, for the configuration of the 
six- vertex model with DWBC which corresponds to A, v{A) is simply the number of vertex 
configurations of type (1), or equally the number of vertex configurations of type (2). 

If A is a permutation matrix, then y{A) is the number of inversions in the permutation tt 
given by (5^.j = Aij. Accordingly, for any ASM A, ^{^) is referred to as the number of 
generalized inversions in A. The statistic viA) was first defined and used by Robbins and 
Rumsey [125, Eq. (18)], who referred to it as the number of positive inversions in A [ , 
p. 182]. A closely-related statistic, Ylii<i<i' <n; i<j'<j<n ^ij ^i'j' = '^{^)+K^)^ was previously 
defined and used by Mills, Robbins and Rumsey [103, p. 344], and provides an alternative 
generalized inversion number for each A G ASM(n). 

Certain generalizations of the statistics (3)-(4) will be defined in Section (3.5). 

It can be checked easily that transposition or anticlockwise quarter-turn rotation of an 
ASM give another ASM, and that, for each A G ASM(n), the statistics (3)-(4) behave under 
these operations according to 

uiA) = uiA^) = !±^- uiA^) - ^(A), /i(A) = /i(A^) = /i(A«), 

Pt(A) = pl(A^) =n-l- pl(A«), pr(A) = pb(A^) = n-l- pt(A«), 

Pb(A) = pr(A^) =n-l- pr(A«), pl(A) = pt(A^) =n-l- pb(A«), (6) 

where Q denotes anticlockwise quarter-turn rotation, i.e., Afj = Aj^n+i-i- By combining 
transposition and anticlockwise quarter-turn rotation, the behavior of any of the eight oper- 
ations of the dihedral group acting on ASMs can be obtained. 

The properties of ASMs with a 1 in a corner can be described easily. For example, for 
the case of ASMs with a 1 in the top-left corner, the sets {A G ASM(n) | An = 1} = {A G 
ASM(n) I Pt{A) = Pl{A) = 0} and ASM(n — 1) are in bijection for any n > 2, where an 
ASM from the first set is mapped to the second set by simply deleting the first row and 
first column. Furthermore, if A from the first set is mapped to A' in the second set, then 
u{A') = z/(A), p(A') = fi{A), pr(A') = pr(A) and pb(A') = pb{A). 

The properties of ASMs in which a 1 on a boundary is separated from a corner by a 
single zero, i.e., ASMs in which a boundary statistic of (4) is 1 or n — 2, can also be 
described relatively easily. For example, for the case of ASMs A with ^421 = 1, the sets 
{A G ASM(n) I A21 = Ai^k+i = 1} = {A e ASM(n) | pl(A) = 1, pt(A) = k} and 
{A G ASM(n -1) \ Au = ... = Ai^k-i = 0} = {A e ASM(n - 1) | pt(A) > A; - 1} are in 
bijection for any n > 2 and k = 1, . . . ,n — 1, where an ASM A from the first set is mapped 
to the second set by replacing A2^k+i by ^2,^+1 + 1, while leaving all other entries unchanged, 
and then deleting the first row and first column. Furthermore, if A from the first set is 
mapped to A' in the second set, then z/(A') = i'{A) — ^2,^+1 — 1 and fi{A') = p{A) + A2^k+i- 
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2.2. Generating functions. In this section, various ASM generating functions are defined, 
and some simple relations which they satisfy are derived. Each of these generating functions 
will be labeled by a certain order of refinement, which corresponds to the number of boundary 
statistics of (4) (or simply the number of boundaries) with which it is associated. Each 
generating function is also associated with the two bulk statistics of (3), i.e., the numbers of 
generalized inversions and — I's in an ASM. 

For each positive integer n, define a quadruply-refined ASM generating function, which 
involves all six statistics of (3)-(4), and associated indeterminates x, y, zi, Z2, z^ and 2:4, as 

^n I3;,y, 2:1, 2:2, 2:3,^4; — 2^AeASM(n) ^ V ^1 ^2 ^3 ^4 • I'J 

It follows that X and y can be regarded as bulk parameters or weights, and that zi, z^^ z^ 
and ^4 can be regarded as boundary parameters or weights. 

Examples of the quadruply-refined ASM generating function, for ra = 1, 2, 3, are 

Zi'"'"^(x, y] zx, Z2, 2:3, za) = 1, 

Z^'^'^'^^x, y; zi, Z2, Z3, Z4) = 1 + x zi z^ 2:3 ^4, 

^quad/ Nil 1 |2 22,222 , 

Z3 (X, y] Zi, Z2, -23, Z4) = 1 + X Zi Z4 + X Z2 Z3 + X Zi Z2 Z^ Z4 + X Z^ Z2 Z3 Z4 + 

3 2 2 2 2 I /o\ 

X z^z^z^z^+xy Zi Z2 Z'i ^4, (8) 

where the terms are written in orders which correspond to those used in (2). 

Now define triply- refined, adjacent-boundary doubly- refined, opposite-boundary doubly- 
refined, singly-refined and unrefined ASM generating functions as, respectively, 

Z*"(x, y; ^1, Z2, zs) = Zr^ix, y- z,, 1, z^, z,) = Ea^asmm ^"^''^ ?/"^^^ 4^^^^ 4"^^^ 4'^^\ 



v(A)„,^l{A) ^Pt(A) PL (A) 



Zf J(x, y- zi, Z2) = Zr'(x, y, z,, 1, 1, Z2) = EAeASM(n) ^'^^^ y'^^^ 4^ ' ^2 

Zr (x, y; zi, Z2) = Zr\x, y- zi, 1, z^, 1) = E^.ASMCn) ^'^^^^ V'^^^ ^f"^^^ 4"^^\ 

Z„(a;, y; z) = Zr'^lx, y- z, 1, 1, 1) = EAeASM(n) ^^^""^ 2/"^^^ ^''"^^^ 

Z^{x,y) = Zr\x,y- 1,1X1) = EAeAsuin)^"^''^ V'^^^ (9) 

where z is a further indeterminate. 

Also define alternative quadruply-refined and alternative adjacent-boundary doubly- refined 
ASM generating functions as, respectively. 



Zr\x,y-z,,Z2,zs,z,) = {Z2z,r~'zr''{x,y;z,,j-,z,,j- 



22 ' -J ' 24 ' 



_ Y^ „u(A) „,p(A) ^Pt(A) n-pR(A)-l pb(^) n-pi^(A)^l 

— l^AeASM(n) X y Zi Z2 Z^ Z4 , 

Zf J(X, 1/; ^1, Z2) = Zr'ix, y; Z„ Z2, 1, 1) = EA.ASM(n) ^^^^^ V'^^^ ^f^^^ 4^^^^- (10) 

Note that Z!^^'^'^{x,y; zi, Z2, z^, z^) is a generating function in which the positions of the I's 
in the first and last columns of an ASM are measured relative to the opposite ends of the 
columns to those used in (4)-(5), i.e., in this generating function, the statistics associated 
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with Z2 and z^ are, respectively, the numbers of O's above the 1 in the right-most column, 
and below the 1 in the left-most column of an ASM. 

It can be seen immediately that some relations among generating functions of (9)-(10), 
involving specializations of boundary parameters to 1, are 

^n"(a;> y\ ^i> i> ^2) = ^^'^■'(a;, y\ zi, Z2), 

Z'/\x, y- 1, zi, Z2) = Zf^ix, y- z^, z^), 
Z*"(x,y;;2i,;22,l) = ^°PP(x,l/;Zi,Z2), 
Zf J(x,|/; z, 1) = Zf J(x,y; z, 1) = Z°PP(x,i/; z, 1) = Z„(x,y; z), 

Zn{x,y]l) = Zn{x,y). (11) 

Some elementary identities satisfied by the ASM generating functions of (7), (9) and (10), 
which follow from the ASM properties outlined in Section 2.1, will now be obtained. 

Note that, in this and all subsequent sections, many of the identities which contain the 
positive integer n will be valid only for all n > 2, or for all n > 3. This will often be due to 
their containing terms (e.g., Zn-i{x,y) or Z„_2(x, ?/)) which are not defined if n is taken to 
be 1 or 2. 

Note also that many of the identities of this section, obtained here using simple combi- 
natorial arguments, can alternatively be obtained as special cases of more general identities 
which will be derived in subsequent sections using other methods. 

Some symmetry relations, which can be derived by acting on ASM(n) with transposition 
or anticlockwise quarter-turn rotation, and using (6), are 

^^""•^(x, y; zi, Z2, Z3, zi) = Z^""'^(x, y; z^, Z3, Z2, zi) 



X 



<-~^y\z,z2z,z,T-'zr'{ii:i.i.i.i) 



Zr\x,y-z,,Z2.z^,z,) = (z,z,z,z,r-^Zr''{x,y;j-^,jr,,j;,j-,) 

Z:\x, y- zi, Z2, z,) = x-(""i)/2 {z,Z2Z,r-' Zf{\, f ; ^, J^, ^; 
Z-f\x,y-z^,Z2) = Z^\x,y-Z2,z^) 



X 



n{n-l)/2 



Z1Z2T ^^n'^Hx'fS^'^)^ 



Zrix,y;z,,Z2) = Z:^^x,y;z2,z,) 



X 



n(n-l)/2 



2i;Z2)""'^°PP(-,^;-,-), 



Z„(x,y)=x"("-i)/2z"'iZ„(i,f). (12) 
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Some identities involving specializations of boundary parameters to 0, which follow from 
the properties of ASMs with a 1 in a corner, as discussed near the end of Section 2.1, are 

ZT'"^ix, y; zi, 0, 2:3, Z3) = Z^^\x, y; zi, 0, ^3) = Zf_^^{x, y; zi, Z3), 

Zf^{x,y;z,0) = Z:P^{x,y;z,0) = Z^_^{x,y;z), 

ZfJ(x,2/;;2,0) = Z„(a;,2/;0) = Z„_i(x,y). (13) 

It will sometimes be useful to refer to boundary parameter coefficients in the adjacent- 
boundary doubly-refined and singly-refined ASM generating functions. In particular, define 

Zf^{x,y)k,M = coefficient of z^'z^^ in Zf^{x,y; Zi, Z2), 

Zn{x, y)k = coefficient of z'^ in Z„(x, y; z). (14) 

These coefficients satisfy identities which correspond to identities of (11)-(13). For exam- 
ple, 

Efca^io ZfK^, y)k,,k2 = Zn{x, y)k„ Efc=J Zn{x, y)k = Zn{x, y) , 
ZfK^,y)k,M = ZfKx,y)k,M, Z4x,y), = x"^"-^)/^ Z„(i, |)„_i_,, 

Zf^ix,y)k,o = Zn^iix,y) 6k,o, Zf^{x,y)k+i,n-i = x'^~^ Z„_i(x,y)fc, 

Zn{x, y)o = Zn-iix, y). (15) 

Furthermore, the properties of ASMs in which a 1 on a boundary is separated from a 
corner by a single zero, as discussed at the end of Section 2.1, lead to identities such as 



Zn^\x,y)k,i = xZn-i{x,y)k^i + y ^^"1=1 Zn~l{x,y)^ - yZn-i{x,y)5k, 



0, 



^n-2 



Zn{x,y)i = xZn-iix,y) +y Yl=i kZn-i{x,y)k. (16) 

3. Previously-known and Related Results 

In this section, an account is given of results for the exact enumeration of n x n ASMs 
(for arbitrary, finite n), involving any of the six statistics of (3)-(4), or any of the associated 
generating functions of (7), (9) or (10). Each of these results is either previously-known, or 
can be obtained straightforwardly from previously-known results. 

Various cases for the values of the parameters x and y, associated with the bulk statis- 
tics (3), are ffist considered separately: y = in Section 3.1, y = x + 1 in Section 3.2, 
X = y = 1 in Section 3.3, and x and y arbitrary in Section 3.4. Further results are then 
discussed in Sections 3.5-3.15. 

3.1. Bulk parameter y = 0. The case in which the bulk parameter y is corresponds 
simply to the enumeration of permutation matrices, with prescribed values of the inversion 
number of associated permutations, and prescribed positions of I's on the boundaries of the 
matrices. 

Let x-numbers and the x-factorial be defined, as usual, as [n]^ = 1 + x -|- . . . + x"~^, and 
[n]J = [n]^[n- 1]^... [1]^. 
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For y = 0, the quadruply- refined ASM generating function is explicitly 

Zr''{x,0;z,,Z2,z,,z,) = 

^0<i<j<n-3^ 



X Z1Z2Z3Z4 l^Q^.^.^^_^[X Z^Z^ +X ' Z^ Z2) 



{xz^zi [n-2]^.24 [n-2]^^^ + ziZ2{xz3Zi)'''^ [n-2]^^^ N-2]xz2 + 

XZ2Z3 [n-2\^^,^ h-2]x23 + Z3Zi{xziZ2Y'^ [n-2]^^,^ [n-2\^^^ [n-?,]^. + 

(1 + x^''~\z,Z2Z3Z,r-^) [n-2]J, (17) 

and tlie ASM generating functions of (9) are explicitly 

Z^\X,Q-Z,,Z2.Z,) = {xZ2Y.,<.^,<n-3{^''^'-'-^AA''''''rX-^^^''zr-^zl) + 

[n-2]^^j + ^2(3^2:1)""^ \n-2\^^^ XZ1Z3 [r2-2]^.^3 [n-?,]J. + 
{zi{xz2Z3T-^ [n-2U, + XZ2 \n-2\^,^ + 1 + ^^^"-^(ziZa^s)""') [n-2]J 
ZfJ(x,0;zi,Z2) = XZXZ2 [n-1]^,, [n-l]^.,^ h-2]J + [n-l]J, 
Zr (x, 0; zi, Z2) = Eo<^<J<n-l (x"+-^--^ ^i z--'-' + x^^^^'^ z^'' zi) [n-2]J, 
Zn{x,0;z) = [n]xz [n-l]xl, 

Z„{x,0) = [n]J. (18) 

Each of these formulae is either a standard result for permutations, or a straightforward 
variation of such a result, and each formula can be derived using simple combinatorial ar- 
guments. See, for example, Stanley [132, Cor. 1.3.13] for a derivation of the last equation 
of (18). 

Note that if (17) is shown first to be valid, then each equation of (18) can subsequently be 
obtained by setting certain boundary parameters in (17) to 1. Alternatively, if the last two 
equations of (18) are shown first to be valid, then (17) and the remaining equations of (18) 
can subsequently be obtained by using the general results (79), (81) (or (84)), (85) and (86), 
which will be given in Section 4. 

In interpreting (17), note also that the seven main terms on the RHS of (17) (i.e., the term 
ending in [n— 4]a,!, the four terms ending in [n — 3]a;! and the two terms ending in [n — 2]^!) cor- 
respond to sums over the sets of n x ra permutation matrices A with {An, Ain, Ann, Ani) equal 
to (0, 0, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1), (1, 0, 0, 0), (0, 1, 0, 0), (1, 0, 1, 0) or (0, 1, 0, 1), respectively. 
In subsequent formulae for quadruply-refined ASM generating functions (e.g., (19), (48) 
and (79)), the RHS will again consist of seven or eight main terms, but these terms no 
longer correspond simply to sums over sets of ASMs with fixed values of the four corner 
entries. 

3.2. Bulk parameters satisfying y = x + 1. The case in which the bulk parameters are 
related by y = x + 1 is closely related to x-determinants (or, in the notation usually used, 
A-determinants) of matrices, domino tilings of an Aztec diamond, tournaments, and the free 
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ferinion case of the six-vertex model. (Note that the sub-case x = 1 and y = 2, which 
corresponds to the so-called 2-enumeration of ASMs, is often considered separately.) Due 
to certain combinatorial and algebraic simplifications, it is again possible to obtain explicit 
formulae. 

For y = X + 1, the quadruply- refined ASM generating function is given explicitly by 

{xZ4Zi + Z4 + Zi — l){xZiZ2—XZi—XZ2—l) X 

{xZ2Zs + Z2 + Zs-l){xZ3Z4-XZ3-XZ4-l) Z^^'^'^{x, X + 1; Zi, Z2, Z3, Z4) = 

;2i22^3^4(x;2i23 + l)(a;22^4 + l)((xzi + l)(xz2 + l)(xz3 + l)(xz4+l))""'(a;+l)("-')("-'^)/'- 

({z2-l){z3-l)z4Zi[{xZ4 + l){xZi + l)y ^ + 

{z4-l){Zi-l)z2Zs{{xZ2 + l){xZ3 + l)y''^j X 
{XZ1Z2 - XZ1-XZ2 - 1){XZ3Z4 - XZ3 - XZi - 1) (x-h l)^""^^^"""^^^^ - 

((;23-l)(^4-l)^1^2(^3^4r-'((x^l + l)(xZ2 + l))""' + 

{Zi-l){z2-l)z3Z4{ZiZ2y''^{{xZ3 + l){xZi+l)y~'^^ X 

(x24^i + Z4 + ^i-l)(a:^2^3 + ^2 + ^3-l)a;"(x+l)("-3)("-^)/2 + 
(2;i-l)(z2-l)(-23-l)(-24-l)((a;-2i-22-a;2:i-a;-22-l)(a;-23-24-a;-23-a^24-l) + 

(xZ4^1 + Z4 + ^l-l)(x^2^3 + ^2 + ;23-l)(^1^2;23^4)""'a;2"-l)(x + l)("-2)("-3)/2, (19) 

and the ASM generating functions of (9) are given explicitly by 

{xZiZ3 + Zi+Z3-l){xZ2Z3-XZ2-XZs-l)Z^^^'{x,X + l; Zi, Z2, Z3) = 

- Z^Z3{xZiZ2 + l){ixZ, + l){xZ2 + l)y~\xZ3 + l)"-\x + lY''-''^^''-^^/^ + 
(22-l)(^3-l)(a;^l23 + ^l + 23-l)^l(^2^3)""'(x2l + l)""V(x + l)("-2)("-3)/2_ 

(^l-l)(23-l)(x^2^3-X^2-X^3-l)(x^2 + l)""'(x+l)("-')("-')/^ 

(xZi^2 + ^l + ^2-l)^fJ(x,X+l;Zi,Z2)=^1^2((x^l + l)(x^2 + l))""'(x+l)("-')("-=^^/' 

- (;^l-l)(^2-l)(x + l)("-l)("-2)/^ 
Z°PP(x,x + l;zi,Z2) = {xziZ2 + l) {{xzi + l){xz2 + l)y'\x + lY^~^^^''''^/\ 

Zn{x,X + l-z) = {xZ+ir'\x + lf'-^^^''-^^'\ 

Z„(x,x+1) = (x+l)"("-^)/^ (20) 

Each of these formulae is either a previously-known result, or a straightforward variation of 
such a result. A derivation of these formulae, based on a result of Robbins and Rumsey [125, 
Sec. 5] for x- determinants, will now be given, since this derivation involves the use of a 
quadratic relation satisfied by the quadruply-refined ASM generating function, and thereby 
shares some features with new material which will be presented in Section 4. An alternative 
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derivation, based on the Izergin-Korepin formula for the partition function of six-vertex 
model with DWBC and the Cauchy double alternant evaluation, will be given in Section 5.8. 
For an n X n matrix of indeterminates, (Ma),^. .^ , define 

Zn{x, M) = Ea.asmh ^"^^^ (^ + 1)^^^^ n:,=i(M.,)^S (21) 

this being the so-called x-determinant of M, as introduced by Robbins and Rumsey [125, 
Sec. 5]. It can be seen that Zn{—1,M) is the standard determinant of M, since the RHS 
of (21) is then a sum over all n x n permutation matrices A, with I'iA) being the number of 
inversions in the permutation associated with A. 

It follows from a result of Robbins and Rumsey [125, Sec. 5] that 

Z„(a;,M)Z„„2(x,Mc) = Z„_i(a;, Mtl) ^„-i(x, Mbr)+xZ„„i(x, Mtr) Z„_i(x, Mbl), (22) 

where Mtl, Mtr, Mbr and Mrl denote the (n — 1) x (n — 1) submatrices corresponding 
to the top- left, top-right, bottom-right and bottom- left corners of M, respectively, and Mc 
denotes the central (n — 2) x (^ — 2) submatrix of M. For x = —1, (22) is the Desnanot-Jacobi 
determinant identity, which will be discussed in more detail in Section 5.5. 
Taking M in (22) to be 



/ 1 z, ... z^^ (^i^2)"-'\ 

M = 



Z4 1 ... 1 zr" 



^n-2 



1 Z2 



(23) 



\izsz,r-' z^-' ... z, 1 / 



it follows, using definitions from (7), (9), (10) and (21), that the quadruply-refined ASM 
generating function satisfies the quadratic relation 

Z^"^'^(x, x+1; zi, Z2, zs, Zi) Zn-2{x, x + 1) = Z'^\{x, x+1; Z4, zi) Z'^\{x, x + 1; 2:2, z-^) + 

xziZ2Z:>,Zi Z^%{x, x + 1; Zi, Z2) Z^%{x, x + 1; Z3, Z4). (24) 

Note that a quadratic relation for arbitrary x and y will be given in (79), and that a different 
quadratic relation for the present case can be obtained by setting y = x + 1 in (79). 
Setting zi = Z2 = Z3 = Z4 = 1 in (24) gives 

Zn{x,X+l) Zn-2{X,X+1) = (x+ 1) Z„_i(x, X + 1)^, (25) 

which, together with Zi{x,x + 1) = 1 and Z2{x,x + 1) = x + 1, gives the formula in (20) for 

Zn{x,X + l). 

Setting ^2 = 2^3 = ^4 = 1 in (24), and relabeling zi as z, gives 

Zn{x,x + 1; z) Zn-2ix,x + l) = (xz + l) Zn^i{x, X + 1; z) Zn-i{x,x + l), (26) 

which, together with Zi{x,x + 1; z) = 1 and the formula for Z„(x,x+1), gives the formula 
in (20) for Zn{x,x+1; z). The formula for Z^^^{x,x + 1; Zi,Z2) can be obtained similarly, by 
setting ^2 = -24 = 1 in (24), and using the formulae for Z„(x, x + 1; z) and Z„(x,x+1). 
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Setting Z2 = Z3 = 1 in (24), and relabeling z^ as Z2, gives 

Zf^{x,x + l;zi,Z2) Z^.2{x,x + 1) = 

Z'^\{x,x + 1; zi, Z2) Zn-i{x, x + 1) + XZ1Z2 Zn-i{x, x + 1; zi) Z„_i(x, x + 1; Z2), (27) 

which, together with Z^ \x,x + l; zi, Z2) = 1, gives 

where, in the sum over i, Zq{x, x+1) is taken to be 1. Substituting the formulae for Zi{x, x+ 
1; z) and Zi{x,x + 1) into (28), and then performing the sum over i and simplifying, gives 
the formula in (20) for Z^'^^{x,x+1; Zi, Z2). 

Finally, the formulae in (19)-(20) for Z^^^'^{x, x+1; zi, Z2, Z3, Z4) and Z^\x, x+1; zi, Z2, z^,) 
can be obtained using (24), together with the formulae for Z^'^^{x, x+1; zi, Z2) and Zn{x, x+1), 
and a similar formula for Z^'^^{x,x + 1; zi,Z2) (which can be obtained from the formula for 
Z^'^^{x,x + 1; Zi, Z2) and the relevant equation from (12)). 

For further information regarding various aspects of the case y = x + 1 and closely re- 
lated cases, including details of several alternative methods for obtaining formulae such 
as (19)-(20), see, for example, Bogoliubov, Pronko and Zvonarev [2U, Sec. 5], Bousquet- 
Melou and Habseiger [21, Sec. 5], Bressoud [23], Brualdi and Kirkland [26], Brubaker, Bump 
and Friedberg [28, 29], Bump, McNamara and Nakasuji [30], Chapman [33], Ciucu [36], [37, 
Thm. 6.1], Colomo and Pronko [45, Sees. 4.2 & 5.2], [46, Sec. 4.3], Elkies, Kuperberg, Larsen 
and Propp [63, 61], Eu and Fu [6ri], Ferrari and Spohn [66], Hamel and King [7~, Cor. 5.1], 
Kuo [91, Sec. 3], Kuperberg ['Jo, Thm. 3], Lascoux [97], McNamara [101], Mills, Robbins 
and Rumsey [103, Sec. 6], Okada [Ul, Thm. 2.4(1), third eq.], Rosengren [126, Sec. 9], 
Striker [131, Sec. 5], [136, Sec. 6], and Tokuyama [^^~, Cor. 3.4]. 

3.3. Bulk parameters x = y = 1. The case in which the bulk parameters x and y are 
both 1 involves numbers of ASMs with prescribed configurations on certain boundaries, but 
without prescribed values of bulk statistics. The cases of zero, one or two boundaries will 
be considered together first, followed by the cases of three or four boundaries. 

Define unrefined, singly-refined, opposite-boundary doubly-refined and adjacent-boundary 
doubly-refined ASM numbers as, respectively. 



An= ASM(n)|, 








An,k= {AeASM{n) 


1 ^i,fc+i = 


= 1}, 




•^ZIm = i^ e ASM(n) 


Ai,ki+i '■ 


^n,n—k2 ~ 


= 1}|, 


<Lm = {^ e ASM(n) 


Ai,ki+i '■ 


= ^fca+l.l = 


= 1}|, 



(29) 
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for < k, ki, k2 < n — 1, with the numbers being for k, ki or k2 outside this range. These 
numbers are therefore related to generating functions of (9)-(10) by 



Zn{l,l)=An, 



zr(i,i;^i,^2) = E:-U-4n 



n-1 yjopp fci fc2 



and to the boundary-parameter coefficients of (14) by 

A,. = ^n(l, 1)., ^:l,,, = ^f ^(1, l)k,M- (31) 

The ASM numbers of (29) satisfy various elementary identities related to their definitions, 
symmetry properties of ASMs or properties of ASMs with a 1 in a corner, and which corre- 
spond to identities, such as (11)-(13) and (15), satisfied by the generating functions. Some 
examples of these identities are 

2—/k=0 '^'n,,k •^ri) •^n,k •^n.n—l—k^ "^71,0 •^n—li 

En— 1 yiopp V-V71— 1 yiadj A 

k2=Q-^n,kiM ~ 2^k2=0'^nMM ~ -^".fci' 

/JOPP _ /JOPP _ /JOPP /jadj _ ^adj 

•^n,kiM •■^n,k2,ki -^n,n-l-ki,n-l-k2^ "^riMM •^n,k2,ki^ 

•^n,fc,0 ~ -An-l Okfl, •^n,kfi ~ "^n.fc+l.n-l ~ -^n-l.fc- (32) 

Furthermore, some identities related to properties of ASMs in which a 1 on a boundary is 
separated from a corner by a single zero, and which can be obtained from (16) and the first 
two identities of (32), are 

yiadj _ spn-2 A , . - A , S, n A , — A 1+ V""^ k A , , — - A 1 

Proceeding to more general identities, for which derivations using combinatorial arguments 
are not currently known, formulae which give the ASM numbers of (29) explicitly are 

•^n — llj=o (n+i)! ' ^ ^ 

r (n+fc-1)! (2n-fc-2)! T-rn-2 (3i+l)! n < U < r, - ^ 

A ^ J fc! (n-fc-1)! (2n-2)! 1 li=0 (n+i-1)!' U \ ft. \ /4 ± 

[0, otherwise, 

jopp _ 1 sr^min(ki,n-k2-l) ( A a , a A — 

•^n,k\,k2 An-i ^i=0 y-^Tijki—i •^n-l,k2+i ' '^n-l,ki-i—l •^n,k2+i 

•^n,ki—i—l'^n—l,k2+i '^n—l,ki~i—l-^n,k2+i+l)y W^/ 

An-l, ki = k2=0 

<Lk2 = { ('t-T') A-1 - E-U E -li i''tT')^^i,n-p l<k,M<n-l (36) 
0, otherwise. 
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It follows from (33)-(34) that the unrefined and singly-refined ASM numbers satisfy simple 
recursion relations, such as 



-'! 



(n-l)!(3n-2)! 



A _ {n-L):(m-z)'. A 
•^" (2n-2)!(2n-l)! -^"-1' 

A A _ 3(3n-2)(3n-4) .2 

-^n-^n-2 4{2n-l)(2n-3) -^n-l' 

k{2n-k-l)An,k = {n-k){n + k-l)An,k-h (37) 

and that the singly-refined ASM generating function at x = y = 1 can be expressed in terms 
of the Gaussian hypergeometric function as 



Zn{l,l;z) = An-12Fl 



1 — n, n 
2-2n ' ^ 



(38) 



and satisfies the hypergeometric differential equation 



z{l-z)-^Zn{l, 1; z) + 2(l-r2-z)^Z„(l, 1; z) + n(n-l)Z„(l, 1; z) = 0. (39) 

The singly- refined ASM numbers also satisfy, for < A; < n — 1, 

A,. = Elo(-l)^rt-7Vn.- (40) 

Recursion relations which involve the opposite-boundary and adjacent-boundary doubly- 
refined ASM numbers are 

( /1°PP _ /1°PP \ A — A A — A A — 

^n—l,fci — l "^71,^2 — 1 T '^n— l,fci — 1 '^n,A;2) v^-'-/ 
•^nM-'^M + '^n,fei,fe2-l ~ -^nMM ~ -^n.fei-l.n-fe ~ ("fci.l ~ l)A:i,o) (0^2,1 ~ L)fe2,o) ^n-1- (42) 

It follows that, in terms of generating functions, (41) can be written as 

{Z1-Z2) ^°P''(1, 1; Zi,Z2) An-l = (zi-l) Z2 Z„(l, 1; Zi) Z„_i(l, 1; Z2) - 

zi {Z2-1) Z„-i(l, 1; zi) Z„(l, 1; Z2), (43) 

and (42) can be written as 

(^1+^2-1) ^:'^J(i,i;^i,^2) = ^i4^r(i,i;^i,i)-(^i-i)(^2-i)A-i. (44) 

The relation (44) will be used in Section 4.3 for the derivation of results for the quadruply- 
refined ASM generating function sX x = y = 1. 

The opposite-boundary doubly-refined ASM generating function aX. x = y = 1 can also be 
expressed as 

,n-l 



Zr (1, 1; ^1, Z2) = 3-"("-i)/^ {q\z, + q){z^ + g))'^-\ 

S(„-l,n-l,...,2,2,l,l) (^+^5 ^f^5 l^^- -jl ) |q=e±2^z/35 (45) 



2n-2 
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where S(„_i^„_i^,..^2,2,i,i)(^^' V+^' ^j ■ ■ ■ ' l) i^ the Schur function indexed by the double- 
staircase partition (n — 1, n — 1, . . . , 2, 2, 1, 1), evaluated at the 2ra parameters ^^' T^' 
1, . . . , 1. Setting zi = Z2 = I in (45), and using (30) and (32), it follows that 

^n = 3^"-("-^^)/2 X (number of semistandard Young tableaux of shape 

{n — l,n — 1, ... ,2,2,1,1) with entries from {1, . . . , 2n}^ . (46) 

The formula (33) can be obtained from (46) using the hook-content formula for semistandard 
Young tableaux. 

A further identity satisfied by the opposite-boundary doubly-refined and unrefined ASM 
numbers is 

n<. t^< MnLkJ = (-l)"("^^)/^+^ (X-l)"-^ (47) 

Proceeding now to the enumeration of ASMs with prescribed configurations on three or 
four boundaries, the alternative quadruply-refined ASM generating function at x = y = 1 is 
given by 

{z4Zi-Z^ + l){ZiZ2-Zi + l){z2Z3-Z2 + l){z3Zi-Z3 + l) Z^'^'"^{1,1; Zi, Z2, Zs, Z^) = 

Zi Z2 Z3 Z4^ deti<ij<4[z^ (zi — l) ~^ Zn-j+i{l,l; Zi)) 

^n—l^n—2^n—3 lll<j<i<4\^* ^j ) 

(z2-l){z3-l)iz4Z,-z, + l){ziZ2-Zi + l){zsZ,-Zs + l)iz2Z4r-'zf\{l,l;j-^,Zi) + 

{Zs~l){z4-l){ZiZ2-Zi + l){z2Z3-Z2 + l){z4Zi~Z4 + l){ZiZsr'^Z';f\{l,l;j-^,Z2) + 
iz4-l)izi-l)iz2Zs-Z2 + l)iz3Z,-Z3 + l)iziZ2-Zi + l)iz2Z4r-'Z^''\il,l;j-^,Zs) + 

(z,-l)iz2-l)izsZ4-Z3 + l){z4Z,-Z4 + l)iz2Zs-Z2 + l)iz,Zsr-'Zl%{l,l;j-^,Z4)- 

{zi-l){z2-l){z3~l){z4-l){{ziZ2-Zi + l){z3Z4-Z3 + l){z2Z4T-^ + 

{z2Z3-Z2 + l)iz4Zi-Z4 + l){ZiZ3r-^)An-2- (48) 

Setting Z2 = 1 in (48) (and then relabeling Z3 as Z2 and z^ as z^), it follows that the 
triply-refined ASM generating function at a; = y = 1 is given by 

{z,Zs-Zs + l){z2Z3-Z2 + l)z^-'Z':\l,l;z^,Z2,j-^) = 

zi Z3 deti<ij<s{z/~'^ {zi-lf~^ Zn-j+i{l, 1; Zi)) 



^n—l^n—2 lll<j 



+ 



<j<3\^i % 



{z2-l){z3-l){ziZs-Z3 + l)zizl^-^Zn-l{l,l;Zi) + 

(z,-l){z3-l)iz2Zs-Z2 + l)z^-'Z^^^,{l,l;Z2). (49) 

Note that the identities (43) and (44), satisfied by the doubly- refined ASM generating 
functions at x = y = 1, can be regarded as special cases of the identity (49) satisfied by the 
triply-refined ASM generating function. More specifically, setting 2:3 = 1 in (49) gives (43), 
while setting Z2 = 1 in (49), and using (43), gives (44). 
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Note also that, by using (44) to replace each case of an adjacent-boundary doubly- refined 
ASM generating function in (48) by an opposite-boundary doubly-refined ASM generating 
function, (48) can be restated as 

2;i 2:2 ^3 ^4 deti<jj<4(Zj {zi — 1) ^^ Zn-j+i{l,l; Zi)) 

'^n — 1 ■^n— 2 -^n— 3 lll<j<7<4V^« ^j ) 
{Z2-1){Z3-1){ZIZ2-ZI + 1){Z3Z4-Z3 + 1)Z4ZIZ2~^Z°^^-^{1,1; Z4, Zi) + 

{Z3-1){Z4-1){Z2Z3- Z2 + 1){Z4ZI- Zi + l)ZiZ2Z^'^Z°^^-^{l,l; Zi, Z2) + 
{Z4-1){ZI-1){Z3Z4-Z3 + 1){ZIZ2-ZI + 1)Z2Z3Z2''^Z°^_^^{1,1;Z2,Z3) + 
{Zi-l){z2-l)iz4Zi- Z4 + 1){Z2Z3- Z2 + 1)Z3Z4Z'^~'^Z°^_^^{1,1; Z3, Z4) + 
{ZI-1){Z2-1){Z3-1){Z4~1){{ZIZ2~ZI + 1){Z3Z4-Z3 + 1){Z2Z4T-^ + 

{Z2Z3 -Z2 + I) {Z4Z1 -Z4 + I) {ziZ3)''-^)An-2. (50) 

The first terms (i.e., the determinant terms) on the RHS of (48)-(50) are proportional 
to certain cases of a function which will be defined in (69). Specifically, using (69), the 
first terms on the RHS of (48) (or (50)) and (49) are ziZ2Z3Z4Xn{l,l; zi,Z2,Z3,Z4) and 
ziZ2Z3Xn{l,l; zi,Z2,Z3), respectively. It will also be seen, in (73), that the function (69) at 
X = y = 1, and hence the first terms on the RHS of (48)-(50), are related to a certain Schur 
function. 

Further expressions for Z^^^'^{1,1] Zi, Z2, Z3, Z4) and Zl^\l,l] Zi,Z2,—), which differ from 
(48)-(49) in the first terms on each RHS, will be obtained in Corollaries 9 and 10. 

The origins of the results given in this section will now be outlined. 

The product formula (33) was conjectured by Mills, Robbins and Rumsey [102, 103, 
Conj. 1], and first proved by Zeilberger [145] and, shortly thereafter, but using a different 
method, by Kuperberg ["^'^j. The product formula (34) was first proved by Zeilberger [116] 
and confirms the validity of further conjectures of Mills, Robbins and Rumsey [102, 103 
Conj. 2]. 

Alternative proofs of (33)-(34) have been given by Colomo and Pronko [15, Sec. 5.3], [ li) 
Sec. 4.2], Fischer [68], and Stroganov [J 10, Sec. 4]. (See also Razumov and Stroganov [117 
Sec. 2], [I 21, Sec. 2] for additional details related to the third of these proofs.) 

The formulae (33)-(34) also follow from certain other known results. For example, as 
already indicated, (33) can be obtained from the result (46) and the hook-content formula 
Alternatively, (33) can be obtained as a special case of a result of Rosengren [126, Cor. 8.4] 
(See also, for example, Rosengren [127] for further related results.) Finally, (34) can be 
obtained by combining a result of Behrend, Di Francesco and Zinn- Justin [17, Thm. 1] 
that An,k is the number of descending plane partitions with each part at most n and exactly k 
parts equal to n, with a result of Mills, Robbins and Rumsey [lU'J, Sec. 5] which gives a 
product formula for the number of such objects. 
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Among these various derivations of (33)-(34), all make essential use of the Izergin-Korepin 
formula [81] for the partition function of the six- vertex model with DWBC, or related prop- 
erties of integrability, except for Zeilberger's proof [145] of (33), which makes essential use 
of a result of Andrews [ ] for the number of totally symmetric self-complementary plane 
partitions in a. 2n x 2n x 2n box, and Fischer's proof [68], which makes essential use of an 
operator formula obtained by Fischer [67, 69] (see also Colomo and Pronko [19, Eq. (3.3)]). 

The hypergeo metric function expression (38) is given by Colomo and Pronko [13, Eq. 
(2.16)], [15, Eq. (5.43)], [16, Eq. (4.19)], and the differential equation (39) is given by 
Stroganov [1 10, Eq. (26)]. 

The relation (40) was obtained by Fischer [68, Sec. 3], with a corresponding relation for 
descending plane partitions, containing a further parameter associated with the sum of the 
parts of a descending plane partition, having been obtained previously by Mills, Robbins 
and Rumsey [102, Sec. 5]. It was also shown by Fischer [6S, Sec. 4] that the formula (34) 
follows from (40), the first three equations of (32), and ^i o = 1- 

The relations (41)-(44) for the opposite-boundary and adjacent-boundary doubly-refined 
ASM numbers were first obtained by Stroganov [140, Sec. 5]. 

The formula (35) follows easily from (41). The formula (36) was first obtained by Fis- 
cher [71, p. 570], and then Ayyer and Romik [ , p. 4], and can be derived by dividing both 
sides of (44) hy zi + Z2 — I and equating coefficients of Zi^Z2' on each side. 

The Schur function expression (45) was obtained by Di Francesco and Zinn- Justin [59, 
Eqs. (2.2) & (2.4)], using a result of Okada [iii, Thm. 2.4(1), second equation]. A general- 
ization of (45) will be given in (73), and a derivation of (73) will be given in Section 5.10. 

The relation (46) was first obtained by Okada [ ' , Thm. 1.2 (Al)]. 

The identity (47) was obtained by Biane, Cantini and Sportiello [18, Thm. 1], by com- 
bining (45) with a certain Schur function identity, also obtained by Biane, Cantini and 
Sportiello [18, Thm. 2]. 

Some further aspects of the opposite-boundary doubly-refined ASM numbers are discussed 
by Fischer [69, Sec. 6]. 

Identities for the quadruply- and triply- refined ASM generating functions at x = y = 1, 
which appear to be equivalent to (48)-(49), were obtained recently by Ayyer and Romik [J 2, 
Thms. 1-2]. The forms of (48)-(49) given here were observed by Colomo [ lU], and can be 
derived by combining a result which will be given in (71), with a result that the partition 
function of the six-vertex model with DWBC is symmetric in all of its spectral parameters 
at its so-called combinatorial point. A derivation of (48) using this approach will be given 
in Section 5.10. (Roughly speaking, this derivation involves associating spectral parame- 
ters ti, ^2, ^3 and ^4 with the first row, last column, last row and first column, respectively, 
of an ASM. The symmetry of the partition function in all of its spectral parameters then 
enables ^2 and ^4 to be associated instead with rows, so that the result (71), which involves 
parameters associated with several rows, can be applied.) 

Note that, since setting certain boundary parameters in (48) to 1 gives (43), (44) and (49), 
Section 5.10 also provides derivations of these other identities. (For the identity (44), such 
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a derivation essentially depends only the symmetry of the partition function in all of its 
spectral parameters, and not on (71).) 

In the versions of (48)-(49) obtained by Ayyer and Romik [12, Thms. 1-2], slightly different 
quadruply-refined and triply-refined ASM generating functions are used, and matrices are 
given which seem to be related by column operations to those in the first terms on the 
RHS of (48)-(49). In particular, the quadruply-refined generating function used there is 

(at X = y = 1), which can be shown straightforwardly to be related to generating functions 
used here by 



Zi Z2 z^ za^ Z^{x, y; zi, Z2, z^, z^) 



ZT'^'^ix, y; zi, Z2, z^, z^ - {z^z^f ' Z^z^^x, y; f-, Zi, 



-1 ^adj 



ZiZ^ ^(XZ3 



f-' Z,f\(x, y; ^1, j-) - z^zr\xzrY-' Z, 
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. n-1i 



-1 r^adj 



{z2Z^Y ' ZlZi{x,y-f- z^] 






((^2^4)"-' + x^^~%z,zsr-') Z„_2(x, y). (51) 



In the matrices given there, functions jniz) and 6n{z) are used which appear to be related 
(although this has not been fully confirmed) to the functions used in the matrices here by 



/ X _ 4(«-l)!(3«-5)! ^2 
/"V^/ (2n-3)! (2n-4)! 



Z ^„_2(1,1;^) (n-3)! (3n-5)! ' (^^1) Znil,l;z), 



X (^\ _ 4(n-l)!(n-4)!(3n-5)!(3n-8)! ^3 y /1 1 . ^N 9(9n^-30n+20} / -, x2 7 H 1 ■ r^ ^^^9^ 



((2n-4)!)2 (2n-5)! (2n-7) 



3.4. Arbitrary bulk parameters x and y. The case in which the bulk parameters x and y 
are arbitrary is of primary interest in this paper, and several new results will be presented 
in Section 4. In this section, the previously-known results for this case, which involve the 
unrefined, singly-refined and opposite-boundary doubly-refined ASM generating functions, 
are reviewed. 

A determinant formula for the opposite-boundary doubly-refined ASM generating function 
is 

(53) 



Z^'^^ix,y;zi,Z2) = det (Knix,y; Zi, Z2)ij), 

0<i,j<n—l 



where 



K„{x,y;zi,Z2] 



V 






irSi'V)-'y 



i—k 



i-l\ (n-l-2\ h^i-k J+1 



2 ' 



^..+1+ ELoEto(::Dr.:7')-V 

^2^k=0 1^1=0 Z^m=Q \i-k) \ k-l H y 



j < n — 3 

J =n-2 (54) 



'-'"z'^z^-'^, j = n-l. 



Setting Z2 = 1 OT zi = Z2 = I'm. (53)-(54), and using standard binomial coefficient identities, 
it follows that determinant formulae for the singly-refined and unrefined ASM generating 
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functions are 






Zn{x,y;z)= det {Kn{x,y;z)i. 

0<«,J<n-l 




Znix,y) = ^^d(^^^^{Kn{x,y)ij), 



(55) 
(56) 

where 

Kn{x, y; z)ij = Kn{x, y; z, l)ij 

K4x, yh = Kn{x, y; 1, 1)., = -5„+i + E^f "^'^ tD et>V-^ (58) 

There are further simphfications ii x = y = 1. For example, 

ir„(l,l),, = -5„+i+(*+^). (59) 

The opposite-boundary doubly-refined ASM generating function satisfies 

(^1-2:2) {zs- Z4) Z°^P{x,y; zi, Z2) Z^P'P{x,y; Z3, Z4) - 

izi-Z3) {z2- Z4) Z°P'^{x,y; zi, Z3) Z^P^^ix^y; Z2, Z4) + 

{z,-z,){z2-Zs)Z^^^{x,y;zuZ^)Z:P^x,y;z2,Z3) = 0, (60) 

and it can be expressed in terms of singly-refined and unrefined ASM generating functions 
as 

{Z1-Z2) Z°^P{x, y; zi, Z2) ^„_i(x, y) = {zi - 1) Z2 ^„(x, y; Zi) Z^-iix, y; Z2) - 

zi {z2-l) Zn-iix,y; zi) Zn{x,y; Z2). (61) 

The two previous identities are essentially equivalent, since (61) can be obtained from (60) 
by setting z^ = 1 and Z4 = 0, and then applying relations from (11) and (13), while (60) 
can be obtained from (61) by expressing each of the six cases of {zi — Zj)Z°^'^{x, y; Zi, Zj) in 
terms of singly-refined and unrefined ASM generating functions, and then checking that the 
resulting expression for the LHS of (60) vanishes. 

The formula (53) was derived by Behrend, Di Francesco and Zinn- Justin [16, Eqs. (21)- 
(22)], as part of the proof of an equality [16, Thm. 1] between the opposite-boundary doubly- 
refined ASM generating function and a four-parameter generating function for descending 
plane partitions. (See also Section 3.12.) The formula (55) was derived by Behrend, Di 
Francesco and Zinn- Justin [17, Eqs. (97)-(98)], prior to the derivation of (53), as part of 
the proof of an equality [17, Thm. 1] (conjectured previously by Mills, Robbins and Rum- 
sey [103, Conj. 3]) between the singly-refined ASM generating function and a three-parameter 
generating function for descending plane partitions. 

For an alternative version of (53), involving a transformation of the matrix Kn{x, y; Zi, Z2), 
see Behrend, Di Francesco and Zinn- Justin [16, Eqs. (65)-(66)], and for alternative ver- 
sions of (55) or (56), involving transformations of the matrices Kn{x,y;z) or Kn{x,y), see 
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Behrend, Di Francesco and Zinn- Justin [17, Props. 1 & 4, Eqs. (65), (66), (87) & (88), & 
Sec. 4.1]. For additional information on determinants closely related to those of (55)-(56), 
including formulae for their evaluation or factorization in certain special cases, see, for exam- 
ple, Andrews [2], Andrews and Stanton [ ], Ciucu, Eisenkolbl, Krattenthaler and Zare [.is, 
Thms. 10-13], Ciucu and Krattenthaler [:')9], Colomo and Pronko [11, Eqs. (23)-(24)], [42, 
Eqs. (4.3)-(4.7)], de Gier [53, Sec. 2.1], Gessel and Xin [7(i, Sec. 5], Krattenthaler [88, e.g., 
Thms. 25-37], [;.w. Sec. 5.5], Lalonde [94, Thm. 3.1], Mills, Robbins and Rumsey [102], [103, 
p. 346], [:i' , Sees. 3-4], Robbins [124, Sec. 2], and Rosengren [128]. 

The case of (61) with x = y = 1 (as given in (43)) was first obtained by Stroganov [1 10, 
Sec. 5], and this relation for arbitrary x and y was first obtained, as an identity involving 
one- and two-point boundary correlation functions for the six-vertex model with DWBC, by 
Colomo and Pronko [11, Eq. (5.32)], [46, Eq. (3.32)]. An alternative proof of (61) was given 
by Behrend, Di Francesco and Zinn- Justin [16, Sec. 5]. This identity will also be obtained 
from Theorem 1 of this paper, as Corollary 4, and it will be found to be a special case of the 
identity (71) given in Section 3.5. 

3.5. ASM enumeration involving statistics associated with several rows or several 
columns. In this section, a generating function involving statistics associated with several 
rows (or several columns) of ASMs is defined, and (in (71)) a result is stated which, for the 
case of particular assignments of certain parameters, provides a determinantal expression for 
this generating function in terms of singly-refined and unrefined ASM generating functions. 
A derivation of this result will be given in Section 5.9. 

For 1 < i,j < n and A G ASM(n), define ASM statistics, associated with row i or column j 
of A, as 



"■^ 



z.--^(A) = \{l<j<n I U^,A,, = E]'=iAAl 
z.™'^(A) = |{1 < z < n I J^,A,, = E]'=iA,,}\, 

f/"'"''{A) = number of — I's in row i of A, 

/i™ •'(A) = number of —I's in column j of A. (62) 

It follows from the defining properties of ASMs that if, for an ASM A, Y17=i^i'j ~ 
Y.f=i^ij' (o^^ equivalently X]i'=i^i'i = Z]j'=i^u')) then Ay = 0. Hence, iy''°'"^{A) and 
i/^^^J(^A) can be regarded as the numbers of certain O's in row i and column j, respectively, 
of A. As will be seen in Section 5.2, for the configuration of the six-vertex model with 
DWBC which corresponds to A, i'^°^''{A) and i>'^°^^{A) are simply the numbers of vertex 
configurations of types (1) and (2) in row i and column j, respectively, of the grid. 

It can also be seen that the statistics z/''°™*(A) and i/™^''(yl) can be written in various other 
ways, for example as 

^coi, ^^) = e:Li (E-;i e^'=i A'^a^ + el. ep=,+i a^.a,). (63) 
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It can be checked easily that, for any A G ASM(ri), the statistics of (3)-(4) and those 
of (62) are related by 



z/"-i(^) = PTiA), z/-'"(A) = pr(A), z/™-"(A) = pb(^), i^^^H^) = Pl{A), 



T.lzir-{A) = EUp-'^'iA) = /i(A), 

^rowl^^) = ;x=°l"(A) = /i™"'^(A) = /i^°'l(A) = 0, (64) 

and that, for any A G ASM(n) and 1 < i < n, the statistics of (62) behave under the 
operations of transposition or anticlockwise quarter-turn rotation of ASMs according to 

zy™"^(A) = i^=°'^(A^) = n - 1 - u^°^\A^) - 2/i''°"*(A), 

fx'^-^'iA) = fx'°^'{A^) = /i"°i'(A«), (65) 

using the same notation as in (6). 

For any < m < n and 1 < ki < . . . < km < n, define an ASM generating function 
associated with rows (or, as will be seen in (68), columns) ki, . . . ,km as 

m 
AeASM(n) i=l 

for indeterminates x, y, zi, . . . , Zm and wi, . . . , Wm- 
It can be seen that 

Z^'--'^"'{x, y;zi,..., Zi_i, 1, Zj+i, . . . , z^; Wi, . . . , Wi_i, y, Wi+i, ..., Wm) = 

^k^,...M-i,k,+i,-,km(^^^ y-zi,..., Zi_u Zi+i, . . . , 2^; u;i, . . . , Wi_i, Wi+i, . . . , w^), (67) 

and that Z^^'---'^'^{x,y;zi, . . . ,Zm',wi, . . . ,Wm) is the unrefined ASM generating function 
Zn{x,y) for m = 0, the singly-refined ASM generating function Zn{x,y; zi) for m = 1 
and fci = 1 or fci = n, and the opposite-boundary doubly-refined ASM generating function 
Z°PP(x, y; Zi, Z2) for m = 2, fci = 1 and k2 = n. 
Using (65), it also follows that 



7 1\/~\ ..... f\/fi 



{x,y]Zi,...,Zm]Wi,...,Wm) 



m 



v(A)„,M(A)~Er=iM^°"=^(^) TT ^^^"'"'(A) ^"'^{A) 



^J- 



AeASM(n) i=l 

^n(n-l)/2( )n-l^fci,...,fc„n^ i ^ . ^ ^. i£^ -™ ) . (gg) 
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Now define, for < m < n, a furtlier function 

ZrXx,y), m = 

Xn{x,y;zi,...,z^) = ^deti<ij<m{z/''^{zi-l)"'-^ Zn-j+i{x,y;zi)) (69) 



1 < m < n. 



Ui<i<j<mi^i-Zj) UT=i Zn-i{x,y) 

Tliis definition is based on a definition of Colomo and Pronko [17, Eq. 6.6], [ , Eq. (4.18)] 
of a closely-related function for the six-vertex model with DWBC 

Some properties of Xn{x, y; zi, . . . , Zm), which follow straightforwardly from (69) and the 
latter equations of (11), (12) and (13), are 



^mj 



Xn{x, y;zi,..., Zm) is symmetric in ^i, . . . , 2;„ 

-^n{X,y,Zi,..., ZjYi) X yZi . . . Zm) "V-r' z' z^' ' ' ' ' ^m / ' 

-^n\Xj y, Zi, . . . , 2^j— 1, i, ■Z^j+l, • • • ; ■2^m/ -^n\X, y, Zx, . . . , 2^j— l, Z^^i, • • • i Zm)i 

Xn\x, y; Zi, . . . , Zj_i, 0, Zi^i, . . . , Zm) = X„_i(a;, y; zi, . . . , Zi_i, 2;j+i, . . . , Zm)- (70) 
The main result of this section is that, for any < m < n and 1 < ki < . . . < km < n, 

^^'■■■'^'"(a;,2/; zi,..., z^; xz^ + {y-x-l)zi + l, . . . ,xZm + {y-x~l)zm+l) = 

Xn{x,y;zi,...,Zm)- (71) 

A derivation of (71), based on a result of Colomo and Pronko [ !~, Eq. (6.8)], [ '", Eq. (A. 13)] 
will be given in Section 5.9. 

Note that, using (67) and the third property of (70), the cases of (71) with m < n all 
follow from the case m = n. 

It can now be seen, using (66) and (71), that the function (69) is a polynomial in x, y and 
zi, . . . ,Zm with integer coefficients. 

It follows from (71) that, for fixed n and m, the LHS of (71) is independent of fci, . . . , km- 
This implies, for example, that the singly-refined and opposite-boundary doubly-refined ASM 
generating functions can be written, for any 1 <k < n and 1 < ki < k2 < n, as 

Znix,y;z) = Z^{x,y;z;xz'^ + {y-x-l)z+l), 
Z°P^{x,y;zuZ2) = Z^^''''{x,y;z^,Z2;xz,^ + {y-x-l)z^ + l,xzi + {y-x-l)z2 + l). (72) 

It can also be seen that (71) contains the identity (61), satisfied by the opposite-boundary 
doubly-refined ASM generating function, as the case m = 2, ki = 1 and ^2 = n. 

It will be shown in Section 5.10 that the function (69) ai x = y = 1 can be expressed as 

X„(l, 1; Zi, . . . , Zm) = 3-'^("-l)/2 (_g)-(«-l) ((^,+g) . . . {Zm + q)^ X 

S(„,-i,„-i,.. .,2,2,1,1) (i7+^) ■■■! i;;rr^' ^2-:_^) lg=e±2-/3' I '<^J 

2n—m 

where this uses the same notation as, and can be seen to generalize, the identity (45). 
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3.6. ASMs with several rows or columns closest to two opposite boundaries pre- 
scribed. In this section, the enumeration of ASMs with prescribed configurations of several 
rows or columns closest to two opposite boundaries (but not involving any further statistics) 
is discussed. 

Let Ki and K2 be subsets of {1, . . . ,n} with \Ki\ + \K2\ < n, and let Ak^ and Ak2 be 
matrices with n columns and \Ki\ or \K2\ rows respectively, in which each entry is 0, 1 or —1, 
along each row and column the nonzero entries alternate in sign, in each row the entries sum 
to 1 , in column j of Aki the entries sum to 1 if j G Ki and the entries sum to with the 
first nonzero entry (if there is one) being 1 if j ^ Ki, and in column j of Ak2 the entries 
sum to 1 if j € K2 and the entries sum to with the first nonzero entry (if there is one) 
being —1 if j ^ K2. (For example, the entries of these matrices could be taken as simply 
{AKi)ij = Sj^Kii and {AK2)ij = ^j,K2i^ where Kii and i^2,j are the ith smallest elements 
of Ki and K2, respectively.) For a fixed choice of matrices A^i and Ak2 which satisfy these 
conditions, define 

•A.n,Ki,K2 = number oi n x n ASMs whose first \Ki\ rows are given by Aki, 

and last |i^2| rows are given by Ak2, (74) 

where it can be checked easily that An^Ki,K2 is independent of the choice of A^i and Ak2- 
(Note that if i^i = or K2 = 0, then there is no restriction on the first or last rows, 
respectively, of the ASMs in (74).) Thus, An,Ki,K2 is the number of n x n ASMs in which 
certain rows or columns closest to two opposite boundaries are prescribed. 

Alternatively, An,Ki,K2 can be written, without reference to A^i and A^^, as the number 
oi {n—\Ki\ — \K2\) X n matrices in which each entry is 0, 1 or —1, along each row and column 
the nonzero entries alternate in sign, in each row the entries sum to 1, and in column j the 
entries sum to 1 if j G {1, . . . , n} \ {Ki U K2), the entries sum to —1 if j G Ki fl K2, the 
entries sum to with the first nonzero entry (if there is one) being —1 if j G Ki\ K2, and 
the entries sum to with the first nonzero entry (if there is one) being 1 if j G K2 \ Ki. 

It can be seen that 

An,Ki,K2 = An,K2,Kx = Ari,{n+l-k\k&Ki} ,{n+l-k\k&K2} (75) 

and that 

'^n,0,0 = An^ An,{k+l},ib = An,k, An,{ki+l},{n-k2} = •^nMM'' ^^^^ 

for any Q < k,ki,k2 < n — 1. 

It can also be seen that ^n,{i,...,n}\Xi,0 or An,{i,...,n}\K2,% give the numbers of possible matri- 
ces Aki or Ak2i respectively, which satisfy the conditions outlined above. Furthermore, for 
a fixed set K of positive integers, ^n,{i,...,n}\x,0 is independent of n, for n > max(_ft'), with 
an operator formula for this number having been obtained by Fischer [uV, Thm. 1]. Some 
related functions for the six-vertex model with DWBC have been studied by Colomo and 
Pronko [49]. 

The numbers An,Ki,K2 have been studied by Fischer [70, 71], Fischer and Romik [72], 
and Karklinsky and Romik [6-->]. Results which have been obtained include linear relations 
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between the numbers An,K,0 with \K\ = 2 and the opposite-boundary doubly-refined ASM 
numbers (Karkhnsky and Romik [83, Eq. (7)], and Fischer [70, Eq. (5.3)]), exphcit formulae 
for An^Kfi with \K\ = 2 (Karklinsky and Romik [83, Thm. 1], and Fischer [70, Eq. (1.3)]), 
and an expression for An,K,{k} with \K\ = 2 in terms of numbers of n x n ASMs with 
prescribed configurations on three boundaries (Fischer [71, Thm. 1]). Also, related numbers 
have been defined (Fischer [70, Eq. (1.2)], [71, Eq. (2.9)]), and have been shown to give 
the numbers An,Ki,K2 ^^ certain cases (Fischer [70, Thm. 1]), and to satisfy certain linear 
relations (Fischer [70, Eq. (5.1) &; p. 253, first equation], [71, Eqs. (1.5) & (1.6)]). 

3.7. The cases Z„(l,3) and Z„(l,3;z). Results involving Z„(l,3) and Zn{l,3]z), some 
of which were previously conjectured by Mills, Robbins and Rumsey [103, Conj. 6-7], [101, 
Conj. 5], and which correspond to the so-called 3-enumeration of ASMs, have been obtained 
by Colomo and Pronko [13], [15, Sees. 4.4 & 5.4], [ IG, Sec. 4.3], Kuperberg [92, Thm. 2], [O.i, 
Thm. 3], Okada [111, Thm. 2.4(1), fourth eq.], and Stroganov [138]. 

3.8. The case Z„(l, 1;— 1). It was observed by Di Francesco [ , Eqs. (2.7)-(2.8)], and 
shown explicitly by Williams [i ii, Thm. 4], that 

Z„(1,1;-1) = V^, (77) 

where Vn is the number of vertically-symmetric n x n ASMs (i.e., ASMs A with Aij = 
Ai^n+i-j)- Note that for n even, it can be seen easily (for example, using the second equation 
of (32)) that both sides of the equation are zero. For n odd, Vn = 111=7 (6^ — 2)!/(?t, -|- 
2i — 1)!, as conjectured by Robbins [123, 124], and first shown by Kuperberg [ ,, Thm. 2, 
second equation]. It is known that Vn is also the number of descending plane partitions 
invariant under a certain operation (see Behrend, Di Francesco and Zinn- Justin [ : , Sec. 4.2 
& Eq. (102)]), and the number of totally symmetric self-complementary plane partitions 
invariant under a certain operation (see Ishikawa [80, Thm. 7.11(i)]), and it has been con- 
jectured by Di Francesco [57, Eq. (2.8)], [58, Eq. (4.5)] that V^ (up to sign) and V^ are 
obtained when a parameter is set to —1 in certain generating functions for totally symmetric 
self-complementary and cyclically symmetric transpose-complementary plane partitions. 

3.9. ASMs with a fixed number of generalized inversions. Expressions for the coef- 
ficients of x^ in Zn{x, 1), i.e., for the number oi n x n ASMs with p generalized inversions, 
can be obtained using a result of Behrend [] 5, Cor. 14], and are given by Behrend, Di 
Francesco and Zinn- Justin [17, Eqs. (24)-(25)]. Some particular cases oi p, for ASMs with 
fixed values of certain further statistics, are also considered by Behrend, Di Francesco and 
Zinn- Justin [ ;, Sec. 7.1], [ \, pp. 337-338]. 

3.10. ASMs with a fixed number of — I's. For the case of Zn{l,y), an expression for 
the coefficient of y"^, i.e., for the number oi n x n ASMs with m —Is, has been obtained 
by Cori, Duchon and Le Gac [50, 99] (and, for m = 1,2, by Aval [ , Prop. 4]), a certain 
factorization property (previously conjectured by Mills, Robbins and Rumsey [103, Conj. 4 & 
Conj. 5], [105, pp. 50 & 54] and Robbins [121, Sec. 2]) has been established by Kuperberg [92, 
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Thin. 3], [9-:i, Thin. 4, first two eqs.], and a certain operator formula containing a parameter 
which corresponds to the number of — I's in an ASM has been obtained by Fischer [(i9, 
Thm. 1]. 

Results involving the coefficient of y^ in Zn{x,y; z), i.e., for a certain generating function 
for n X n ASMs with a single —1, have been obtained by Lalonde [94, 96]. 

3.11. Objects in simple bijection with ASMs. Simple bijections are known between 
ASMs and various other combinatorial objects. Some examples, all of which are reviewed by 
Propp [113, Sees. 2-4 & 7], are bijections between ASM(n) and sets of (n + 1) x (n + 1) corner- 
sum matrices, (n + 1) x (n + 1) height-function matrices, 3-colorings of an (ri + 1) x (n + 1) grid 
with certain boundary conditions, monotone (or Gog) triangles with bottom row 1, . . . ,n, 
configurations of the six-vertex model on an n x n grid with DWBC, and fully packed loop 
configurations on an n x n grid with certain boundary conditions. Some further examples 
are a bijection between ASM(n) and the set of sets of n osculating paths on an n x n grid 
in which all paths start and end on two adjacent boundaries (see, for example, Behrend [15, 
Sees. 2-4], Bousquet-Melou and Habseiger [21, pp. 68-69], Bressoud [22, pp. 226-227], or 
Egecioglu, Redmond and Ryavec [62, pp. 35-36]), and a bijection between ASM(n) and the 
set of alternating paths for any fixed fully packed loop configuration on a.n n x n grid (see 
Ng [110, Prop. 3.1]). 

In each of these cases, the statistics for the other combinatorial object which correspond, 
under the bijection, to the ASM statistics (3)-(4) or (62) can be obtained straightforwardly. 
For example, this will be done for configurations of the six- vertex model with DWBC in (111). 
Accordingly, all of the results reviewed or obtained in this paper for ASM(n) with the statis- 
tics (3)-(4) or (62) can alternatively be expressed in terms of any of these other combinatorial 
objects and their associated statistics. 

In the case of the six- vertex model with DWBC, ASM generating functions are closely re- 
lated to certain partition functions, expectation values, probabilities or correlation functions 
for the model, and the results of this paper could alternatively be expressed in terms of the 
latter quantities. For example, the quadruply-refined ASM generating function is closely as- 
sociated with a certain four-point boundary correlation function for the six- vertex model with 
DWBC. For studies of various correlation functions, and related quantities, for the six- vertex 
model with DWBC, see, for example, Bogoliubov, Kitaev and Zvonarev [19], Bogoliubov, 
Pronko and Zvonarev [20], Colomo and Pronko [44, 46, 47, 49], Foda and Preston [73], and 
Motegi [106, 107]. 

3.12. Descending plane partitions. Descending plane partitions are certain combinato- 
rial objects first defined by Andrews [-,•>]• It was shown by Behrend, Di Francesco and 
Zinn- Justin [16, Thm. 1] that the opposite-boundary doubly-refined ASM generating func- 
tion Z°PP(x, y; zi, Z2) is equal to the generating function for descending plane partitions with 
largest part at most n, in which the statistics associated with x, y, Zi and Z2 are, respectively, 
the number of nonspecial parts, the number of special parts (as first defined by Mills, Robbins 
and Rumsey [103, p. 344]), the number of n's, and the number of (n — l)'s plus the number 
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of rows of length n — 1, in a descending plane partition. (Note that no further statistics for 
descending plane partitions are currently known which, together with the previous four sta- 
tistics, lead to generating functions which are equal to the adjacent-boundary doubly-refined, 
triply-refined or quadruply-refined ASM generating functions, and that no general bijection 
is currently known between ASM(n) and the set of descending plane partitions with largest 
part at most n.) 

It follows that certain known results for descending plane partitions correspond to results 
involving the opposite-boundary doubly-refined, singly-refined or unrefined ASM generating 
functions. For example, certain such results were outlined in Section 3.4. 

Results are also known for the enumeration of descending plane partitions with a prescribed 
sum of parts. In particular, it follows from work of Mills, Robbins and Rumsey [ , Sec. 5] 
that, for < k < n ~ 1, 

V r,\D\ - rr'^-l [3»+ll.;' 

Z^DeDPP(n) y — llj=0 [n+j]q! ' 

Y^ JD\ _ gfe" [n+fc-1]^! [2n-k-2]gl T-rn-2 [3i+l]gl (^o\ 

Z^DeDPP(n,fc) y [fc]q! [n-fc-1],! [2n-2]q! 1 lj=0 [n+i-1],! ' y'°) 

where DPP(n) denotes the set of descending plane partitions with largest part at most n, 
DPP(?7,, k) denotes the set of elements of DPP(n) with exactly k parts equal to n, \D\ denotes 
the sum of parts of a descending plane partition D, and [n]g\ is the usual g-factorial of n. As 
expected, each RHS of (78) reduces to each RHS of (33)-(34) for q = I. However, no ASM 
statistic is currently known which corresponds to the sum of parts of a descending plane 
partition, and leads to generating functions matching those of (78). 

For further information, and related results or conjectures, regarding descending plane 
partitions, see, for example, Andrews [2, 3], Ayyer [1(J], Behrend, Di Francesco and Zinn- 
Justin [16, 17], Bresssoud [ J], Bressoud and Propp [25], Krattenthaler [90], Lalonde [94, 95, 
96], Mills, Robbins and Rumsey [102, 103, 105], Robbins [123, 124], and Striker [135]. 

3.13. Totally symmetric self-complementary plane partitions. Totally symmetric 
self-complementary plane partitions are certain combinatorial objects first defined by Mills, 
Robbins and Rumsey [iUl, Sec. 1], and Stanley [131, Sec. 2]. It was shown by Fonseca and 
Zinn- Justin [75], following conjectures of Mills, Robbins and Rumsey [ , Conj. 2 & 3], 
that the opposite-boundary doubly-refined ASM generating function with both bulk pa- 
rameters set to 1, Z°^^{1,1] zi,Z2), is equal to a generating function for totally symmetric 
self-complementary plane partitions in a 2n x 2n x 2n box in which the statistics associated 
with zi and Z2 can be certain pairs from among several statistics defined by Doran [til. 
Sec. 7], and Mills, Robbins and Rumsey [104, Sec. 3]. The singly- refined case of this result 
was obtained earlier by Razumov, Stroganov and Zinn- Justin [122, Sec. 5.5]. (Note that 
no statistics for totally symmetric self-complementary plane partitions are currently known 
which have the same enumerative behavior as the number of generalized inversions or the 
number of — I's in an ASM, and that no general bijection is currently known between ASM(n) 
and the set of totally symmetric self-complementary plane partitions in a 2?2 x 2ri x 2n box.) 
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It follows that certain known results for totally symmetric self-complementary plane par- 
titions correspond to results involving the ASM generating functions Z°^^{1,1; Zi,Z2) or 
Zn{l,l;z), or the ASM numbers A^^^,^^, An,k or An- For example, for Z°pp(1, 1; zi, ^2), 
Zn{l, 1; z) or An, Pfafiian expressions follow from results of Ishikawa [79, Thms. 1.2 & 1.4, 
& Sec. 7] and Stembridge [^■'■', Thm. 8.3], constant-term expressions follow from results of 
Ishikawa [79, Sec. 8], Krattenthaler [nV, Thm.] and Zeilberger [145, Sublems. 1.1 & 1.2], and 
integral expressions follow from results of Fonseca and Zinn- Justin [75, Eqs. (4.9) & (4.14)]. 
Several of these results include a certain ASM statistic, first introduced by Mills, Robbins 
and Rumsey [104, p. 287], but not considered in this paper. 

For the case of adjacent-boundary doubly-refined ASM enumeration, it has been conjec- 
tured by Cheballah [ ; , Conj. 4.3.1] that A^l^j^^ equals the number of totally symmetric 
self- complementary plane partitions in a.2nx2nx 2n box for which a certain pair of statistics 
have values fci and ^2- 

For further information, and related results or conjectures, regarding totally symmetric 
self-complementary plane partitions, see, for example, Andrews [1], Ayyer, Cori and Gouyou- 
Beauchamps [11], Bressoud [22, Sec. 6.2], Bressoud and Propp [25], Cheballah [3 1], Cheballah 
and Biane [35], Di Francesco [57, 58], Doran [61], Fonseca [71, Sees. 3-4], Fonseca and 
Zinn- Justin [7 ], Ishikawa [79, 80], Krattenthaler [S7], Mills, Robbins and Rumsey [10 1], 
Robbins [123], Stanley [131], Stembridge [133, Sec. 8], Striker [134, 136], Zeilberger [145], 
Zinn- Justin [1 19], and Zinn- Justin and Di Francesco [150]. 

3.14. Loop models. Numbers of certain ASMs, or plane partitions, have been found to 
appear also as appropriately normalized entries (or sums of such entries) of particular eigen- 
vectors associated with certain cases of integrable loop models (or associated quantum spin 
chains). A wide variety of such cases are known, and these often involve ASMs with pre- 
scribed link patterns of associated fully packed loop configurations, with prescribed values of 
some of the statistics of (3)-(4) or of related statistics or parameters, or subject to invariance 
under certain symmetry operations. 

Of the many papers which study the confirmed or conjectured appearances of numbers of 
ASMs in such contexts, a few examples are Batchelor, de Gier and Nienhuis [13], Cantini and 
Sportiello [31, 32], Di Francesco [55, 56], Di Francesco and Zinn- Justin [ ' ], Pasquier [112], 
Razumov and Stroganov [115, 116, 118], and Stroganov [137]. For reviews of some of these 
matters, see, for example, de Gier [51, 52, 53], or Zinn- Justin [149]. 

3.15. ASMs invariant under symmetry operations. This paper is focussed on the enu- 
meration of ASMs with prescribed values of the bulk statistics (3) and boundary statis- 
tics (4), but with no other conditions applied. However, other studies have focussed on the 
enumeration of ASMs which are invariant under certain symmetry operations, or subject 
to related conditions, and which in some cases also have prescribed values of some of the 
statistics of (3)-(4). For results of this type (some of which remain conjectural), see, for ex- 
ample. Aval and Duchon [8, 9], Bousquet-Melou and Habseiger [21], Bressoud [22, Sec. 6.1], 
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Kuperberg [ ], Okada [111], Robbins [123, 124], Razumov and Stroganov [117, 119, 120], 
Stanley [ < ], and Stroganov [139, 141]. 

Numbers of certain such ASMs have also been related (again, in some cases, only con- 
jecturally) to numbers of certain descending plane partitions (see, for example, de Gier, 
Pyatov and Zinn- Justin [51, Prop. 3, first equation], and Mills, Robbins and Rumsey [ ■ , 
Conj. 3S]), numbers of certain totally-symmetric self- complementary plane partitions (see, 
for example, Ishikawa [79, 80], and Mills, Robbins and Rumsey [101, Conjs. 4 & 6]), or 
appropriately normalized entries of eigenvectors associated with certain cases of loop models 
(see, for example, the references given in Section 3.14). 

4. New Results 

In this section, the main result of this paper is stated, and some of its implications are 
derived and discussed. The proof of the result will be deferred to Section 5.6. 

4.1. Main result. The primary result of this paper is as follows. 

Theorem 1. The quadruply-refined ASM generating function satisfies 

y{zi- Z2){zi- z^) Zl"'"^{x, y; zi, Z2, zs, z^) Zn-2{x, y) = 

{{zi~l){z2-l)-yziZ2){{z3-l){zi~l)-yz:iZ^) Z^\{x,y] z^, zi) Z^\{x,y] Z2, Z'i) - 

{x{zi-l){zi-l)-y){x{z2-l){z^-l)-y) ziZ2Z'iZiZl\{x,y] zi, Z2)Zl\{x,y] Z'i, Zi) - 

{z2-l){z-i-l)[{zi-l){zi-l)-yZiZi) Zfj^{x,y;z4,zi) Zn^2{x,y) + 

{z-i-l){zi-l){x{zi-l){z2-l)-y) ziZ2 {xz'iZiY'^ Z^\{x,y]Zi,Z2) Zn-2{x,y) - 

{zi-l){zi-l){{z2-l){z^-l)-yz2Z^) Zf_l^{x,y;z2,Z3) Z„_2(a;,y) + 
{zi-l){z2-l){x{z3-l){z4^-l)-y) zsZ4,{xziZ2)'^~'^ Z'^\{x,y;z3,Z4) Zn^2{x,y) + 

(^1-1)(^2-1)(^3-1)(^4-1)(1 - ix^ZiZ2ZsZ^r-') Zn-2{x,yf. (79) 

The identity (79) holds for each n > 3. Furthermore, if ZQ{x,y) is taken to be 1, then 
it can be seen, using Z^ \x,y; zi, Z2) = Z^ \x,y; zi, Z2) = 1 and Z2^^ {x,y; zi, Z2, z^, Z4) = 
1 + xziZ2Z^Zi (from (8)-(10)), that (79) also holds for n = 2. 

The proof of (79), which will be given in Section 5.6, will involve using a relation between 
a certain generalized ASM generating function and the partition function of the six-vertex 
model with DWBC, and then combining the Desnanot-Jacobi determinant identity with the 
Izergin-Korepin formula for this partition function. 

An alternative form of (79) will be obtained, and some further related results will be 
discussed, in Section 5.7. 

It can be seen that (79) enables Z'jl^^'^{x,y; Zi, Z2, -23, -24) to be obtained recursively, using 
the initial conditions (from (8)) Z'^^^ {x,y; zi, Z2^ Zs^ Z4) = 1 and Zg"*^ {x,y; Zi, Z2^ z^^ z^) = 
1 + xziZ2Z^Zi, and the definitions (from (9)-(10)) Z^^{x^y;zi,Z2) = Z'jl^^'^{x,y; Zi,l,l,Z2), 
Zfi{x,y;z,,Z2) = ZrH^^V, ^u ^2,1,1) and Z„(x,y) = Zr''{x,y;l,lA,l). Accordingly, 
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in this sense, Z!jl^'^'^{x,y; Zi, Z2, Z3, z^), and all of the ASM generating functions of (9)-(10) 
which are defined in terms of Z^^^'^{x, y; zi, Z2, z^, Z4), are determined by (79). 

Note, however, that if the generating functions are computed recursively in this way, then, 
for each successive n, Z^^^'^{x, y; zi, Z2, ^3, Z4^) should first be computed for arbitrary zi, Z2, ^3, 
and Z4, with the factor {zi — Zs){z4 — Z2) being explicitly canceled from both sides of (79), so 
that division by zero is avoided when boundary parameters need to be set to 1 in subsequent 
computations. Alternatively, expressions in which this cancelation has effectively been done, 
will be given in Section 4.2. 

By replacing Z2 and z^ by — and — respectively. Theorem 1 can be restated for the 
alternative quadruply- refined ASM generating function of (10) as 

y{zi- z^){z2- Zi) Zl"'"^{x,y; zi, Z2, z^, Zi) Zn-2{x,y) = 
{{z,~l){z2-l)+yz,){{z3-l){zi-l) + yz3) (^2^4)""' x 

z^M^, y, ^, ^1) K-ii^, y; ^, ^3) - 

{x{zi-l){zi-l)+yz4){x{z2-l){zs-l)+yz2) 21^3(^2^4)""^ x 

K-i{x, y, zi, ^) K-i{x, y; ^3, ^) - 

{z2-l){z,-l){{z4~l){z^-l)+yz,){z2Z4)^~'z::'\{x,y;j-^,z,)Zr,^2ix,y) + 
{z3-l){zi-l){x{zi-l){z2-l) + yz2) ziz^-^ixzs)"-'^ Zf\{x,y;zi,j^) Zn-2ix,y) - 

{Zi-l){z,~l){{z2~l){z,-l)+yz,) {z2Zi)''-' Z:f\{x,y;j-^,Z3) Zr,.2{x,y) + 

{zi-l){z2-l){x{z3-l){zi-l)+yzi) Z3Z^''^{xzi)''~^ Z'^%{x,y;z3,J2) Zn~2ix,y) + 

(;2i-l)(^2-l)(^3-l)(^4-l)((^2^4)"-' - {x^z.z^)^-') Z^_2{x,yf. (80) 

4.2. Consequences for arbitrary bulk parameters x and y. In this section, some conse- 
quences of Theorem 1 , for the case in which the bulk parameters x and y remain arbitrary, are 
derived. These results are obtained from (79) using only definitions and elementary proper- 
ties of ASM generating functions from Section 2.2, and the initial conditions (from (8)-(10)) 
Zf\x,y;zi,Z2) = Zf\x,y;zi,Z2) = 1 and Zf\x,y; zi, Z2) = Zf\x,y; zi, Z2) = 1 + XZ1Z2. 
In particular, most of the results are obtained by setting boundary parameters to 1, and 
using (11) for the corresponding specializations. 

Alternative forms of some of the results of this section will be given in Section 5.7. 

Corollary 1. The triply-refined ASM generating function satisfies 

{z2-zi){z3-l) Z^J\x,y; zi, Z2, Z3) Zn^2ix,y) = 

{{z2-l){z-i-l)-yz2Z3,) zi Z^%{x,y;zi,Z3) Zn-i{x,y; Z2) - 

{x{zi-l){z3-l)-y) Z1Z2Z3 Zl\{x,y; Z2, Z3) Zn-i{x,y; zi) - 

{zi-l){z-i-l) Z2Zn^i{x,y;z2) Zn-2ix,y) + 

{Z2-1){Z3-1) Zi {XZ2Z3)'^'^ Zn^l{x, y, Zi) Zn-2{x,y). (81) 
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Proof. Set Z2 = 1 in (79) (and then relabel z^ as Z2 and z^ as 23). D 

Corollary 2. The adjacent- boundary doubly-refined ASM generating functions satisfy the 
recursion relations 

(zi- l)(z2-l) Zf\x, y; zi, Z2) ^„_2(a;, y) = yziZ2 Z^\{x, y; Zi, Z2) Zn-i{x, y) + 

[x{Zi-l){z2-l)-y) ZiZ2Zn-l{x,y]Zi) Zn^i{x,y]Z2) + 

{zi-l){z2-l) Zn-i{x,y) Zn-2{x,y), (82) 

{zi-l){z2-l) Z^\x,y] zi, Z2) Zn-2{,x,y) = yziZ2Zf_!^{x,y; Zi, Z2) Zn-iix,y) + 
{{zi-l){z2-l)-yziZ2) Zn^i{x,y; zi) Zn-i{x,y; Z2) + 

(zi-l)(z2-l) {xziZ2)"-'^ Zn-iix,y) Zn-2{x,y). (83) 

Proof. To obtain (82), set Z2 = 1 in (81) (and then relabel z^ as ^2)- To obtain (83), set 
2:1 = 1 in (81) (and then relabel Z2 as zi and z^ as 2:2). □ 

Note that, if Zo{x,y) is taken to be 1, then (82) and (83) hold for all n>2. 
Corollary 3. The triply-refined ASM generating function also satisfies 

y{z2-zi)zs Z^^\x, y; Zi, ^2, Z3) Z„_i(x, y) = 

{zi-l)({z2-l){z3~l)-yz2Z3) Zf^{x,y; Zi, Z3) Zn-iix,y; Z2) - 

{z2-l){x{zi-l){z3-l)-y) ziZsZf^{x,y;z2,Z3) Zn-iix,y;zi) - 

{zi-l){z2-l){z^-l) Zn-i{x,y\ Z2) Zn-i{x,y) + 

{zi-l){z2-l){z^-l) ziz^-\xz:,Y Zn-i{x,y-zi) Zn-i{x,y). (84) 

Proof. Use (82)-(83) to replace Z^_\(x, y; 2:1, z^) and Z^_\(x, y; ^2, ^3) in (81) by terms which 
instead contain Z^^{x, y; Zi, z^) and Z^'^^{x, y; Z2, Z3), and then cancel an overall factor which 
contains a term 2:3 — 1. D 

Corollary 4. The opposite-boundary doubly-refined ASM generating function satisfies 

{Z1-Z2) Z°^'^{x, y; zi, Z2) Zn-\[x, y) = {zi - 1) 2:2 Z„(x, y; Zi) ^„_i(x, y; Z2) - 

zi (2:2-1) Zn-i{x, y; zi) Z„(x, y; 2:2). (85) 

Proof Set Z3 = 1 in (84). D 

Note that Corollary 4 is a previously-known result, as already given in (61) and discussed 
in Section 3.4. 
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Corollary 5. The adjacent-boundary doubly-refined ASM generating functions can be ex- 
pressed as 

n-l 






-V(^l-1)(^2-1) 

{x{zx-\){z2-\)-y) Zi{x,y\z^ Zi{x,y\Z'2) 



yZi_i{x,y)Zi{x,y) 



(86) 



-V(^i-1)(^2-1), 

i-w ^n-i , {zxZ2T~'^^{iz^-^){z2-\)-yzxZ2)Zi{x,y;zi)Z,{x,y;z2) \\ 
" ^"^'^^ + yZ.-i(x,,)Z.(x,,) JJ'('') 

where, in the sums over i, Z^^x^y) is taken to be 1. 

Proof. It can be checked straightforwardly (again taking ZQ{x,y) = 1), that the initial con- 
ditions Zl ■'(x, y; zi, Z2) = Zl \x, y; Zi, Z2) = 1 and recursion relations (82)-(83) are satisfied 
by (86)-(87). D 

Corollary 6. The two types of adjacent-boundary doubly-refined ASM generating function 
are related by 

{{zi-l){z2-l)-yziZ2) Zf^{x,y; Zl, Z2) - {x{zi-l){z2-l)-y) ziZ2Zf^{x,y; zi, Z2) 

= (zi-l)(^2-l)(l - ixziZ2r) Zn-i{x,y). (88) 

Proof. This can be obtained directly from (86)-(87). Alternatively, it can be obtained 
from (82)-(83), by showing that {{zi-l){z2-l) - yziZ2) Z^^\x,y]Zi,Z2) - {zi-l){z2- 
l)Zn-i{x,y) and {x{zi-l){z2-l)~y) Z1Z2 Z'^^^{x,y; zi, Z2) - {zi-l){z2-l){xziZ2Y Zn-i{x,y) 
satisfy the same recursion relation and initial condition. D 

Note that, from (12), the two types of adjacent-boundary doubly-refined ASM generating 
function are also related by 

Zf^{x,y- zi,Z2) = x<--^^l^ {Z1Z2T-' Z^^^l, |; ^, ^). (89) 

Combining (88) and (89), it follows that the adjacent-boundary doubly-refined ASM gener- 
ating function satisfies 

{{zi-l){z2-l)-yziZ2) Z^^\x,y-zi,Z2) = 

{x{z,~l){z2-l)-y) {Z1Z2T x<--'^" Zf J(i, I- ^, ^) + 

{zi-l){z2-l){l-{xziZ2T)Zn.i{x,y). (90) 

By using (79) and (88), together with identities from (12), it is also possible to ob- 
tain a certain identity which involves Z^^'^{x, y; Zi, Z2, Z3, z^) and Z^^'^{x, y; Z3, Z2, Zi, z^) (or 
Z^^^'^{x,y; Zl, Z2, Z3, Z4) and Z^^'^{x,y;zi,Z4^,z^,Z2)). A form of this identity will be given 
in (141). 
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Corollary 7. The boundary parameter coefficients in the singly-refined ASM generating 
function, as defined in (14), satisfy 



Zn(,x,y)k = Zn-i{x,y)6kfi + Zn^i{x,y)^{y''^^i . j I -.^ 



+ 



k—i—l n—i—2 

Zn-i-i{x,y)Zn-i-2{x,y) PV ^-1 )\ ^ 



i v^ v^ Zn-i^i{x,y)j^ Zn-i-i{x,y)j., l(k-ii-2\ fn-j2-2 



jl=0 J2=0 



w*-rO("r'iii'(-) 



where Zo{x,y), if it appears, is taken to be 1. 



Proof. Expanding the factors l/{{zi — l)(-22 — l))" * in (86) as binomial series, equating 
coefficients of z^^Z2^ on both sides of (86), and using the definitions (14), gives 

min(fci,fc2)-l / ,, _ s ., _ X 



i=0 
fci— j— 1 k2—i—l 



i sr^ sr^ Zn-i^i{x,y)j^ Z^_i_i{x,y)j^ ( Ai-ji-2\ /A;2-j2-2 
^0 i^o Zn-i-i{x,y)Zn-i-2{x,y) y\ i-l )\ i-1 



h-n-i\fk2-j2-i\ , , (32) 



Summing (92) over k2, using a standard binomial coefficient summation identity, and rela- 
beling fci as k, then gives (91). D 

Corollary 8. The unrefined ASM generating function satisfies 
Z^{x,y) = Z^^^{x,y)U +J2U^'(^^^lJ + 



(93) 



Zn-i-i{x, y) Z„_i_2(a;, y) 
where ZQ{x,y) is taken to be 1. 
Proof. Sum (91) over fc, and use a standard binomial coefficient summation identity. D 

Note that various other expressions for Zn{x,y)k and Zn{x,y), which provide alternatives 
to (91) and (93), can be obtained by taking the hmits Zi — )■ 1 and 2:2 — )• 1 of (86) in ways 
somewhat different from those used in Corollaries 7 and 8. 

As observed in Section 4.1, the main result (79) enables all of the ASM generating functions 
of (7), (9) and (10) to be computed recursively. The corollaries of this section essentially 
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comprise special cases of (79) which apply to specific ASM generating functions, and can 
alternatively be used for their computation. 

For example, it can be seen that (91) and (93) give Zn{x,y)k and Zn{x,y) in terms of 
Zi{x,y)k and Zi{x,y) for i = 1, . . . ,n — 1, and thereby enable the singly-refined and unre- 
fined ASM generating functions to be computed recursively. The adjacent-boundary doubly- 
refined, opposite-boundary doubly-refined, triply-refined and quadruply-refined ASM gener- 
ating functions can then be computed using (86)-(87), (85), (81) or (84), and (79), respec- 
tively. 

Note that the opposite-boundary doubly-refined, singly-refined and unrefined ASM gener- 
ating functions can instead be computed using the determinant formulae (53), (55) and (56). 

4.3. Consequences for bulk parameters x = y = 1. In this section, some consequences 
of Theorem 1, for the case in which the bulk parameters x and y are both set to 1, are derived. 
In contrast to the consequences of Theorem 1 given in Section 4.2 for the case of arbitrary x 
and y, which were obtained essentially using only (79), the consequences in this section are 
obtained using (79) together with the nontrivial relation (44) between the adjacent-boundary 
and opposite-boundary doubly-refined ASM generating functions at x = y = 1. 

Corollary 9. The alternative quadruply-refined ASM generating function at x = y = 1 
satisfies 

(z4Zi-Z^ + l){ZiZ2-Zi + l){z2Zs-Z2 + l){z3Z4-Zs + l) Z^'^'^'^il,!; Zi, Z2, Zs, Z^) = 
Zi Z2 ^3 ^4 



X 



An-2{ZI-Z^){z2-Zi) 

[{ziZ2-Zi+l){ziZ2-Z2+l){z^Zi-z^+l){z-iZi-Zi+l) Z°^]'^{1, 1] Z4, Zi) Z''^\{1, 1; Z2, 2:3) - 

{ZAZi-Z4+l){ZiZi-Zi + l){z2Zz-Z2 + l){z2Z^-Z-i + l) Z"^\{1, I] Zi, Z2) ^^!.1(1, 1; 2^3, Z^)) + 
{z2-l){z-i-l){ziZ2- ZI + 1){Z2,Z4:- Z^ + 1)Z4:ZIZ2'^ Z'l^\{l,l] Z4:, Zi) + 
{Z'i-l){Zi-l){z2Z'i-Z2 + l){z4^Zi-Zi + l)ZiZ2Z'^~^Z°^]'^{l,l]Zi,Z2) + 
{Zi-l){zi-l){z^Zi- Z^ + l){ziZ2- ZI + 1)Z2Z^Z2''^ Z°^]'^{l,l] Z2, Z-i) + 
{Zi-l){z2-l){ZiZi-Zi + l){z2Z'i-Z2 + l)Z'iZ4^z'l~^Z''^]'^{l,l;Z'i,Zi) + 
{zi-l){z2-l){z^-l){Zi-l){{ziZ2-Zi + l){z^Zi-Z^ + l){z2ZiY-^ + 

{z2Z^-Z2 + l){z^Zi-Zi + l){ZiZ^Y-^)An-2- (94) 

Proof. Setting x = y = I'm. (80), and using (89), gives 

{Zi-Zz){z2-Zi)Z^^'"^{l, 1; Zi, Z2, Z^, Z4) An^2 = 

{z,Z2-Z2 + l){z^z,-z, + l){z2Z,r-^ zf\{\, 1; ^, zy) Zt\{\, 1; ^, z^) - 
(^4^i-^i + l)(z2^3-^3 + l)(^i^3)"-' Z:l\(l, 1; J^, Z2) Zf\{\, 1; ^, z,) - 

(^2-l)(^3-l)(^4^1-^4 + l)(^2^4)""'^.f\(l,l;^,2l)X-2 + 
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(z3-l)(^4-l)(^1^2-^l + l)(^1^3)"-'^:-\(l,l;i^,^2)A-2- 

(^l-l)(^2-l)(^3^4-^3 + l)(^1^3r"'^:'\(l,l;^,^4)A-2 + 

(Zl-1)(^2-1)(^3-1)(^4-1)((^2^4)"-' - {z.Zsr-') Al_,. (95) 

Multiplying both sides of (95) by {z4Zi-Z4+l){ziZ2-Zi+l){z2Z3-Z2+l){z3Z4-zs+l) / {{zi-zs){z2- 
2:4) A-2), and using (44), in the form {ziZ2-Zi+l)z'l~'^Zf^^{l, 1; ^, Z2) = Z2Z°^^^{1, 1; Zi, Z2) + 
z'^-^{zi-l){z2-l)An-2, then gives (94). D 

Corollary 10. The triply-refined ASM generating function at x = y = 1 satisfies 



iz^z,-z, + l)iz2Z,-Z2 + l)z^-'Z':\l,l;z,,Z2,j-J = 

— ^^- — {{ziZ3-Zi + l){ziZ3-Z3 + l)z2 ^„_i(l, 1; zi) Z°P'i(l, 1; ^2, Zs) - 

vA„-2 (^1-^2) (^3-1) 

{z2Zs-Z2 + l){z2Z3-Zs + l)zi Zn-iil,l; Z2) Z°^_^^{1,1; Zi, Z3)) + 

{z2-l){z3~l){ziZ3-Z3 + l)ZizJ^-^Zn-l{lA;Zi) + 

{z,-l){z3-l)iz2Z3-Z2 + l)z^-'Z^^_,{l,l;Z2). (96) 

Proof. Set 2:2 = 1 in (94) (and then relabel Z3 as Z2 and z^ as Z3). D 

It can be seen that (94) and (96) provide alternative formulae to (48)-(50) for the quad- 
ruply- and triply-refined ASM generating functions at x = y = 1. In fact, (94) and (96) 
differ from (50) and (49), respectively, only in the first terms on each RHS. It can also be 
seen that the last five terms on the RHS of (94) could be replaced by the last five terms on 
the RHS of (48). 

5. Proofs 

In this section, proofs of Theorem 1, and of several of the results of Section 3 are given. 
Preliminary results which are needed for these proofs are obtained or stated in Sections 5.1- 
5.5, while the main steps of the proofs are given in Sections 5.6 and 5.8-5.10. In Section 5.7, 
alternative forms of certain results of Section 4, and some further results, are discussed. 

5.1. A generalized ASM generating function. In this section, a generating function 
which generalizes the quadruply-refined ASM generating function is introduced, and some 
of its elementary properties are identified. This generating function will be shown, in Sec- 
tion 5.3, to be proportional to a certain case of the partition function of the six- vertex model 
with DWBC. 

The generalized ASM generating function involves the six statistics of (3)-(4), together 
with four further statistics associated with the entries in the corners of an ASM, and is 
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defined, for indeterminates x, y, Zi, Z2, z^,, z^, Z41, Z12, Z23 and Z34, as 
^T^i^^ 2/5 ^i5 ^2, Z3, Z4] Z41, Z12, Z23, Z34) = 

E^u{A) „ h(A) ^Pt(A) PbXA) pb(A) pl(A) i_Aii 1-Ai„ 1-A„„ ^ 1-A„i /q^n 
•^ y ^1 ^2 ^3 ^4 ^41 ^12 ^23 ^34 • W ' / 

AeASM(n) 

For example, for n = 3, the function is 

-^3 (3^) y; ^1) ^2, 2:3, Zi, Z41, Zi2, Z23, Z34) = Z\2 Z34 + X Z1Z4 Z41 Z\2 Z34 + 

X Z2 Zs Z\2 ^23 ^34 + 2^ Z\ Z2 ^3 ^4 ^41 ^12 ^23 + 2^ ^i ^2 ^3 ^4 ^41 ^23 ^34 + 

X z^ Z2 z^ z^ Z41 Z23 + xy Zi Z2 Z3 Z4 Z41 Z12 Z2Z Z34, (98) 

where the terms are written in an order which corresponds to that used in (2). 

It can be seen immediately that the quadruply-refined ASM generating function is 

Z^'^ix, y, zu Z2, Z3, Z4) = Zf\x, y- z,, Z2, z^, Z4, 1, 1, 1, 1). (99) 

By acting on ASM(n) with transposition or anticlockwise quarter-turn rotation, and us- 
ing (6), it follows that 

^n'^iXj y] Zi, Z2, Z3, Z4] Z41, Z12, Z23, Z34) = Z^^^{x, y; Z4, Z3, Z2, z\] Z41, Z34, Z23, Z12) 

= x"('^-l)/2 {z,Z2Z3Z4r-' Zr(i, f ; ^, ^, ^, j-^; ^12, ^23, ^34, Z41) . (100) 

The properties of ASMs with a 1 in the top-left corner, as outlined in Section 2.1, imply 
that 

■^T'^yXi y'l Zlt ^2; ^3; Z4, 0, 2:12, 2:23, 2:34) = Z^™(X, y] 0, 2:2, 23, Z4, Z41, Z12, Z23, Z34) 

= Z12 Z34 ^f!"i(x, y; 1, Z2, Z3, 1; 1, 1, ^23, 1), (101) 

and the fact that an ASM cannot contain a 1 in both its top-left and top-right corner, implies 
that 

Z^'^ix, y- 2:1, Z2, 2:3, 2:4; 0, 0, 2:23, 2:34) = 0. (102) 

It can be seen from the definition (97) that Z^'^{x, y; zi, Z2, Z3, Z4; Z41, zu, Z23, Z34) is linear 
in 2:41, 2:12, 2:23 and 2:34. Therefore, using interpolation for these parameters at and 1, it 
follows that 

^n^^yXj y] Zl, Z2, Z3, Z4; Z41, Zi2, Z23, Z34) = 

J2 (1 - ^AlY-'^ ^,1' (1 - ^uY""' Z,^- (1 - ^23)'-^^ ^23^^ (1 - ^34)'-^^ Z,t X 

n,i2,i3,i4e{0,i} ..... , s 

Z° {x,y;zi,Z2,Z3,Z4;ii,i2,i3,t3)- (103) 

Applying cases of (9)-(10) and (99)-(102) to (103) (as a result of which, for example, (100) 
and (102) imply that nine of the sixteen terms on the RHS of (103) vanish, while (9) and (99)- 
(101) give Z^^^{x, y; zi, Z2, Z3, Z4; 0, 1, 0, 1) = Zn-2{x, y)), it now follows that the generalized 
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ASM generating function can be expressed entirely in terms of quadruply- refined, adjacent- 
boundary doubly-refined and unrefined ASM generating functions as 

^T'^i^^ V'l ^li ^2; Z^i Z4', Z4i,Zi2, Z23, 234) = 2^41 ^12 2:23 Z34 Z^^^ [X, y, Zi, Z2, Z^, Z4) + 

Zii Zi2 (1-2:23) zzi Zl\{x, y; z^, zi) + 
Z41 Zi2 Z2'i (1-2:34) 2:1^2 {xz'iz^'^'^ ^^-i(a;, y\ Zx, Z2) + 

(1-Z4i) 2:12 2:23 Z34 ^^^{x, y; Z2, Zs) + 
Zai (1-2:12) 2:23 2:34 2:3^4 (x2;iZ2)''"^ K%i^^ y^ ^3, 2:4) + 

(1-2:41)^12 (1-^23) 2:34 ^n-2(a;,y) + 

2:41 (1-2:12) 2:23 (1-2:34) x^"-^ (2:1 2:2 2:32:4)""^ Zn-2ix, y). (104) 
Some special cases of (104), which will be used in Section 5.6, are 

^f "(a;, y; zi, l, l, Z2; z^, l, l, l) = Zf^'^ix, y, 1, zi, Z2, 1; 1, 1, 2:12, 1) 

= Zi2Z^'^^{x,y;zi,Z2) + (1-2:12) Z„_i(x,|/), 
^f "(a;, y, zi, Z2, 1, 1; 1, 2:12, 1, 1) = ^^"(0;, y; 1, 1, 2:1, 2:2; 1, 1, 1, 2:12) 

= 2:12 ZfJ(x,y; 2:1, 2:2) + (1-^12) (x2;i2;2)"-^Z„_i(x,y). (105) 

5.2. The bijection between ASMs and configurations of the six- vertex model with 

DWBC. In this section, the set of configurations of the six- vertex model on an n x n grid 
with DWBC is described, and the details of a natural bijection between this set and ASM(n) 
are summarized. This is standard material, with similar accounts having been given, for 
example, by Behrend, Di Francesco and Zinn- Justin [16, Sec. 5.1], [17, Sees. 2.1 &; 3.1]. 

The six-vertex, or square ice, model is a much-studied integrable statistical mechanical 
model (see, for example, Baxter [11, Chaps. 8 & 9] for further information and references), 
with DWBC for the model having been introduced and first studied by Korepin [85]. The 
bijection between ASMs and configurations of the model with DWBC was first discussed 
by Elkies, Kuperberg, Larsen and Propp [61, Sec. 7], with the details having mostly been 
observed previously, but using different terminology, by Robbins and Rumsey [125, pp. 179- 
180]. There exist closely- related bijections between ASMs and certain sets of osculating 
lattice paths (see Section 3.11 for references), and between ASMs and certain fully-packed 
loop configurations (see, for example, Propp [113, Sec. 7]). 

Let Qn be the nxn undirected grid with vertex set {{i,j) \ i,j = 0, . . . , ?t.-|-1}\{(0, 0), (0,'n,-|- 
1), [n + 1, 0), [n + l,n + 1)}, where {i,j) is taken to be in the ith row from the top and jth 
column from the left, and for which there are horizontal edges between {i,j) and {i,j ± 1), 
and vertical edges between {i,j) and (i ± l,j), for each i,j = 1, . . .n. This grid is shown 
in Figure 1. Each vertex or edge of Qn can be described as either internal or external, i.e., 
the n^ vertices of degree 4 are internal, the remaining 4n vertices of degree 1 are external, 
the 2n(n — 1) edges which connect two internal vertices are internal, and the remaining 4n 
edges are external. 
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(0,1) ■ ■ ■ (0,n) 



(1,1) 



(»,1) 



(l,n) 



(n,n) 



1(1, n+l) 



>(n,n + l) 



(n + 1,1) ..." (n + l,n) 

Figure 1. The grid Gr, 



A configuration of tfie six-vertex model on Qn witli DWBC is an assignment of arrows to 
tfie edges of Qn such that the arrows on the external edges on the upper, right, lower and left 
boundaries of Qn are all directed upward, leftward, downward and rightward, respectively, 
while the arrows on the four edges incident to any internal vertex satisfy the condition that 
two point towards and two point away from the vertex. 

Now define 6VDW(n) to be the set of all configurations of the six-vertex model on Qn with 
DWBC. For example. 



6VDW(3) 




(106) 



In an element of 6VDW(r2), there are six possible configurations of arrows on the four 
edges incident to an internal vertex of Qn- These so-called vertex configurations, are shown 
in Figure 2, where the numbers 1, . . . , 6 will be used to label the types, as indicated. 



>-^'-►- 









>~^^->- 



->~^^^ 



V 

-^-■^-k- 



(1) (2) (3) (4) (5) (6) 

Figure 2. The possible arrow configurations on edges incident to an internal vertex. 



For C G 6VDW(n), let Cij G {1,...,6} denote the type of vertex configuration in C 
at internal vertex (z,j) of Gn, and let ^f(k)iC), ^fl°^'{C), ^f^|^^\C) and A/'|)(C) denote 



m 



the numbers of type-fc vertex configurations in C in the whole of Qn, in row i of Qn, 
column j of Qn, or at vertex (z, j) of Qn, respectively. Thus, for example, AC^^JC) = Scij,k 
and Af^k){C) = |{(z, j) I 1 < ^,J < n, Q, = k}\ = El,=i^U^)- 
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It can be shown easily that, for any C G 6VDW(n), 

Mii){C) = M(2){C), Ar(3)(C) = A/'(4)(C), Ar(5)(C) = M(,){C) + n, (107) 

fl^\C) = Ar(-;"(C) = ^l^riC) = A/'(5) 



Ar;r 1(C) = AfttiC) = Afir'iC) = Af[^\C) = 1. (108) 



(See, for example, Bressoud [22, p. 228], [23, p. 290] for a discussion of the first two equations 
of (107).) It can also be seen that, for any C G 6VDW(n) and I < i,j < n, 

Cy G {1, 3, 5}, Cn G {2, 3, 5}, C^, G {2, 4, 5}, Q, G {1, 4, 5}, (109) 

and therefore that 

CnG{l,5}, Ci„G{3,5}, ane{2,5}, C„iG{4,5}. (110) 

It can be shown straightforwardly that there is a natural bijection between ASM(n) and 
6VDW(n), and that, for each A G ASM(n) and C G 6VDW(?2) which correspond under this 
bijection, the statistics (3)-(4) and (62) satisfy 

iy{A) = Ar(i)(C) (= Ar(2)(C)), f,{A) = Ar^e){C) (= ^^,){C) - n), 
PTiA)=Afl^;\C), pR(A)=Ar(-;"(C), pB(A)=Ar(™)""(C), pUA)=Af[^;\C), 

/i™"'(A)=Ar('g7*(C), fi'°^^{A)=Af^^^'{C). (Ill) 

The details of this bijection are as follows. To map A G ASM(n) to C G 6VDW(?2), for 
each i = 1, . . . ,n and j = 0, . . . , n, assign a right or left arrow to the horizontal edge of Qn 
between {i,j) and {i,j + 1), according to whether the partial row sum Yl''j'=iAi3' i^ o'^ 1; 
respectively. Similarly, for each i = 0, . . . , n and j = 1, . . . , n, assign an upward or downward 
arrow to the vertical edge of ^„ between (i, j) and (i + 1, j), according to whether the partial 
column sum Yll'=i ^i'j is or 1, respectively. (The fact that each partial row sum and each 
column sum is or 1 is given by (1).) To map C G 6VDW(n) to A G ASM(n), let 

1, Cjj = 5 

Aj = { -1, a, = 6 (112) 

0, a, G {1,2,3,4}. 

The behavior of the statistics fi, p^, pr, Pb and pr under this bijection can be seen 
easily. To obtain the behavior of the statistics z/, z/™^* and z/^°'-', note that, for each A G 
ASM(n) and C G 6VDW(ra) which correspond under the bijection, and each 1 < i,j < n, 
X]i'=\A'j = Y.]'=iAij' = (or equivalently Yll'=iAi'j = Zlj'=\^ii' = 0) if and only if Cj^ = 1, 
^^^ Z]i'=\^i'i = Y^f=iAij' = 1 (o^^ equivalently J2l'=iAi'j = I]j'=\^ii' = 1) if and only if 

(-^ij = 2. 

As examples of this bijection, in (2) and (106) the elements of ASM(3) and 6VDW(3) are 
listed in an order for which respective elements correspond under the bijection. 
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5.3. The partition function of the six-vertex model with DWBC. In this section, 
the partition function of the six- vertex model witli DWBC is introduced. A relation between 
this partition function, for certain assignments of its parameters, and the generalized ASM 
generating function (97), for certain assignments of its parameters, is then derived using the 
bijection of Section 5.2. 

Let a weight W(^k){u,v) be associated with the vertex configuration of type k, where u 
and V are so-called spectral parameters. 

The partition function for the case of the six-vertex model of relevance here depends on 
these weights, and on spectral parameters Ui and Vj associated with row i and column j 
of Qn, for each 1 < i,j < n. Specifically, this partition function is defined as 



Ce6VDW(n) i,j=l 

Let the weights now satisfy 

W(^i){u,v) = W(2)iu,v) = a{u,v), W(^s){u,v) = W(^4){u,v) = b{u,v), 

W(^5)iu,v) = W^e)iu,v) = c{u,v), (114) 

for functions a, b and c. 

Taking the spectral parameters in (113) to be 

U2 = ... =Un-l =r, V2 = ... = Vn-l = S, Ml=tl, i^n = ^2, Un = h, Vi = t^, (115) 

for indeterminates r, s and ti, . . . ,^4, and using (107)-(110), gives 



(7e6VDW(n) 

\ a(r,s) / \ b{r,s) / V c{r,s) ) 

('a{r^\^(t "(^)-^r2"(C) (b±i2)\^^l^^^{C)-^fl^){C) /£(r^N l-Ar,V5(C)-Ar-(C) 
V a(r,s) J V b{r,s) / \ c(r,s) / 

\ a{r,s) I \ b{r,s) ) \ c{r,s) ) 

t a(rM) \^m ' (^)-^a) iC) I' b(r^N ^(T; ' {C)-N^^^ (C) / c(rAO \ 1-A^(V) (C)-^("4 (C) 
V a{ryS) ) \ b(r,s) ) V c{r,s) ) 

l a(hM) \^fn,iC) r c(H,U) \J^U^i.C) r b{t„t^) xNl^,{.C) ^ c{t,M) \^I^^C) ^ 
\ a{r,s) ) \ c{r,s) / \ b{r,s) J \ c(r,s) J 

/ a(t3.f2) \Ar-(C) / c(t3,f2) \^r5?(^) ( b(t3M) \-^m(C) (sitsMl\^(f)i(^) 
V a{r,s) ) V c{r,s) ) \ b{r,s) J \ c(r,s) J 



X 



X 
X 
X 
X 



(116) 
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It can be seen that 

Ar(i)(C) + Ar(3)(C) + Ar(6)(C) = ^^, (117) 

Afl^j^^lO + Afi^^'^ic) = Ar[^l\c) + Ar[:l\c) = n-i, (iis) 



Afp,^{c)+Arjiuc) = m]{c)+ArrMc) = i, (119) 



where (117) follows from (107) and ELiA/'(fc)(C) = n^, (118) follows from (108)-(109) and 
ELi -^^"^(C) = ELiA^(T)'(C) = n, and (119) follows from (110). 

By using (117)-(119) to eliminate Ar^3){C), Afl°^\C), A/'(™/"(C), Afl^^^'iC), Af;^°J\C), 
Afl^){C), A/'(^5')(C), Ar(^^(C) and A/'g(C) from (116), and then using the bijection of Section 5.2 
between ASM(n) and 6VDW(?t,), the behavior (111) of the statistics (3)-(4) under this 
bijection, and the definition (97) of the generalized ASM generating function, it now follows 
that (116) can be written as 

Z^^{ti,r,...,r,ts;U,s,...,s,t2) = 

b(r s)"-("--l) c(r s)" ( '>(ti''')'>(r,t2)b(t3,s)b{r,t4) Y-^ c{tuU)c{tut2)c{Ui,t2)c{t3,t4) 



b(r,s)* / c{ti,s)c{r,t2)c{t3,s)c(r,t4) 



X 



rygen I ( a{r,s)\'i / c{r,s)\'2 _ a{ti,s)b{r,s) a{r,t2) b{r,s) a{t3,s)b{r,s) a{r,t4) b{r,s) _ a(r,s) a(ti, 44) c(ti,s) c(r,t4) 
"■ U b{r,s) ) ' \b{r,s)) ' a(r,s) 6(ti,s) ' a(r,s) 6(r,t2) ' a(r,s) fe(t3,s) ' a(r,s) fe(r,t4) ' a(ti ,s) a(r,t4) c(r,s) c(ii ,44) ' 



6(r,8) b(ti ,t2) c(ti ,s) c(r,t2) a(r,3) a(t3,fa) 0(^3,5) c(r,t2) b(r,s) bjtsfy) c{t3,s) c(r,t4) 
b{ti,s) b{r,t2) c{r,s) c{ti,t2) ' a{t3,s) a(r,t2) c{r,s) c{t3,t2) ' b{t3,s) b(r,t4) c{r,s) c{t3,t4) 



. (120) 



5.4. The Izergin— Korepin determinant formula. In this section, the Izergin-Korepin 
formula for the partition function (113), with certain assignments of the weights (114), is 
stated. 

It was shown by Izergin [81], using certain results of Korepin [85], that if the weights (114) 
satisfy the Yang-Baxter equation (see, for example, Baxter [11, pp. 187-189]), then the 
partition function (113) is given by an explicit formula involving the determinant of a certain 
n X n matrix. 

Let the functions a, b and c in (114) be given by 

a{u, v) = u q^'"^ — V q^^'^, 

h{u, v) = V q^''^ — u q'^'"^ = a{v, u), 

c{u,v) = {q-q-^)u'^^v^/\ (121) 

where g is a further indeterminate, often known as the crossing parameter. (Note that q will 
be present in many subsequent equations, even though it may not appear explicitly.) 
It can be seen that these functions satisfy 



a{u, vy + h{u, vY — c{u, v) 
a{u,v) b{u,v) 



2 



■(g + g-i). (122) 
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The fact that the LHS of (122) (which is often denoted as 2A) is independent of u and v 
iniphes that the six-vertex model weights given by (114) and (121) satisfy the Yang-Baxter 
equation (see, for example, Baxter [14, Eq. (9.6.14)]). 

The resulting Izergin-Korepin determinant formula is then given by the following result. 

Theorem (Izergin). The partition function (113), with weights given by (114) and (121), 
satisfies 



Ul<i<j<ni'^i - Uj){Vj - Vi) l<^,J<n \a{Ui,Vj) b{Ui,Vj)J 

This theorem can be proved by showing that each side of (123) satisfies, and is uniquely de- 
termined by, four particular properties. Specifically, each side is symmetric in ui, ... ,Un and 
in f 1, . . . , f„ (which can be obtained for the LHS using the Yang-Baxter equation, and is im- 
mediate for the RHS), is a polynomial of degree 2n — 1 in each of u^ , . . . , Un , f i , . . . , fn , 
satisfies the same recursion relation for cases in which Ui = q^^Vj for some i and j, and sat- 
isfies the same initial condition at n = 1. For further details of this proof, see, for example, 
Izergin [81], Izergin, Coker and Korepin [82, Sec. 5], Fonseca and Zinn- Justin [75, Sec. B.l], 
Kuperberg [92, Sec. 1], or Zinn- Justin [149, Sees. 2.5.2-2.5.3]. For an alternative proof, see 
Bogoliubov, Pronko and Zvonarev [20, Sec. 4]. 

Note that, although the determinant and the denominator of the prefactor on the RHS 
of (123) both vanish if Ui = Uj or Vi = Vj for some i ^ j, the RHS has a well-defined limit in 
these cases, as a polynomial in Ui , . . . , Un , f i , . . . , f n . Accordingly, it will be valid to 
use (123) to derive properties of the partition function for the assignments (115), as will be 
done in Section 5.6. 

5.5. The Desnanot— Jacobi identity. In this section, the Desnanot-Jacobi determinant 
identity is stated, and discussed briefly. 

For a square matrix M, let Mtl, M^r, Mbr, Mbl and Mq be the submatrices of M 
defined in Section 3.2. The Desnanot-Jacobi identity, in the form which will be used here, 
can be stated as follows. 

Theorem (Desnanot-Jacobi). For any square matrix M , 

det M det Mq = det Mtl det Mbr - det Mtr det Mrl- (124) 

For an algebraic proof of (124), see, for example, Bressoud [ , Sec. 3.5], and for a combi- 
natorial proof, see Zeilberger [1 17]. Cases of (124) for n x n matrices M with small values 
of n were pubhshed by Desnanot in 1819 (see Muir [108, Eqs. (A)-(G), (A')-(G'), pp. 139- 
142]). The further attribution to Jacobi is based on the fact that, for arbitrary n, (124) 
corresponds to the case m = 2 of the identity, published by Jacobi in 1834 (see Muir [108, 
Eq. (XX. 4), p. 208]), that for any nxn matrix M and any m < n, each mxm minor of the 
matrix of (n — 1) x (n — 1) minors of M equals the complementary minor of M multiplied by 
(det M)™"^. For proofs of the Jacobi identity, using a variety of methods, see, for example. 
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Brualdi and Schneider [27, Sec. 4], Knuth [m, Eq. (3.16)], Leclerc [100, Sec. 3.2], Muir [109, 
Sec. 175], or TurnbuU [143, pp. 77-79]. 

For a discussion of alternative forms of (124), and of further determinant identities of 
which (124) is a special case, see, for example, Behrend, Di Francesco and Zinn- Justin [l(i. 
Sec. 4]. 

5.6. Proof of Theorem 1. In this section, the main steps in the proof of Theorem 1 
are given. These involve using the relation (120) between the generalized ASM generating 
function and the partition function of the six-vertex model with DWBC, the Izergin-Korepin 
formula (123), and the Desnanot-Jacobi identity (124). 

By applying the Desnanot-Jacobi identity (124) to the matrix ( c^u^v^} — ^ \ ^^^ 

then applying the Izergin-Korepin formula (123) to each of the six determinants which 
appear, it follows that the partition function (113), with weights given by (114) and (121), 
satisfies 

(Mi - Un) (fn - ^^l) Zf^{Ul ■ ■ ■ , Un, Vi, . . . , V^) Z'^_^{u2, . . . , M„_i; t^2, • • • , Vn-l) = 

a{Ui, Vn) b{Ui, Vn) a{Un, Vi) b{Un, Vi) Z^_^{ui, ..., Un-l] fl, . . . , Vn-l) X 

Z^Yi(m2,. . .,Un,V2,. . . ,Vn) 
- a{Ui, Vi) b{Ui, f i) a{Un, Vn) b{Un, Vn) ^^Yi(Mi, • • • , Un-l, ^2, • • • , fn) X 

Znlliu2, ...,Un;Vi,..., Vn-l), (125) 

for any Ui, . . . , u„, Vi, . . . ,Vn- In fact, it can be seen, from the previous derivation, that (125) 
is satisfied by any function which has the form of (123), for arbitrary functions a, b and c. 

Certain forms of the Desnanot-Jacobi identity have previously been combined with cer- 
tain cases of the Izergin-Korepin formula by Korepin and Zinn- Justin [86, Sec. 3] (see also 
Sogo [129, Sec. 4]), and by Behrend, Di Francesco and Zinn- Justin [16, Sec. 5.4]. 

Now let parameters x, y, zi, Z2, z^, z^, z^i, zu, Z23, Z34, z'^^, z[2, z'^j, and ^34 be given in 
terms of parameters g, r, s, ti, ^2, ^3 and ^4 by 

■^ \b{r,.s)J ' y \b(r,s)J ' 

_ a{ti,s)b{r,s) _ a{r,t2) b{r,s) _ a{t3,s)b{r,s) _ a(r,t4,) b{r,s) 

^1 a(r,s) 6(ti,s)' ^2 a{r,s) b(r,t2) ' ^ a(r,s) b(t3,s) ' ^^ a(r,s) 6(r,i4) ' 

_ a(r,s)a(ti,f4) _ b{r,s)b{ti,t2) _ a{r,s) ajtsfy) _ b{r,s)b{t3,t4) 

^41 a(ti,s)a(r,t4)' ^^ b(ti,s)b{r,t2)' 23 a(t3,s) a{r,t2) ' ^4 b{t3 ,s) b{r,U) ^ 

I _ b{r,s)b{ti,t4) I _ a{r,s)a{ti,t2) I _ b{r ,s) bjtz ,12) I _ a{r,s)a(ti,ti) (10P,\ 

^41 — b(ti,s)b(r,t4)' ^12 ~ a(ti,s)a(r,t2)' ^23 " 6(43,5) b(r,t2) ' ^34 " a{t3,s) a(r,t4) ' \^^^) 

where the functions a, b and c are given by (121). 

The parameterizations (126) allow arbitrary x, y and ^i, . . . , Z4 to be expressed in terms 
of g, r, s and ti, . . . ,^4. In particular, it follows, using (121), (122) and (126), that q can 
first be taken to be a solution of 



X 



^/2g2 + (a;-l/ + l)g + x^/2 = 0, (127) 
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that r and s can then be taken to satisfy 

(x^/2 + g)r = (x^/2g + l)s, (128) 

and that ti, . . . , ^4 are then given by 

, _ s(xV^ zi g+1) , _ rjx'^/^ Z2+q) , _ s{x^/^ z^q+l) , _ r(xV^ 24+g) H 90") 

It can be checked easily that the parameterizations (126) imply that Z41, zu, Z23, -234, -241, 
z'12, -223 ^^^ ^u '^^^ t)e expressed in terms of x, y and 2:1, . . . , 2:4 as 

^ _ I iz4-i){zi-i) ^ _ ]^ x{zi-i){z2-i) ^ _ -^ iz2-i)iz3-i) ^ _ ]^ x(za-i]{z4-i) n3Q^ 

41 2^242:1 ' 12 y ) 23 1/22^3 ' 2/ ' V / 

and 

^/ _ j^ x{z4.-l)izi-i) ^ _ -^ (2:i-l)(2:2-l) ^ _ ^ x{z2-V){z-a-\) ^ _ -^ (Z3-1)(Z4-1) ^3^ 
41 y ? 12 y zi 22 ' 23 1^ '34 y 23 24 ' *^ ' 



_ (2i-l)(2j-l) ^^ ^i^y J_ J_ 

The parameterizations (126) also give 



Thus, each Zi^ and z-,- has the form /(y, 2;^, z^) = 1 - ^""^ ^ j^^^^. — - or /(^, ^^ , ^ 



(t2 - U) iU - h) = "^^'^^''^^;,g,y;g,ffi'^^^^'''*^^ {Z4 - Z2) {Zl - Zs). (132) 

Assigning the parameters in (125) according to (115), applying the relation (120) to each 
of the six cases of the partition function in (125), and using (126) and (132), it now follows 
that 

iz4-Z2) {Z1-Z3) Zf''{x, y; zi, Z2, Z3, z^; ^41, Z12, Z23, Z3i) Zf\{x, y\ 1, 1, 1, 1; 1, 1, 1, 1) = 

y Zl Z2 Z3 z^ {zyi z'12 ^34 ^4 ^f^ii^^ Vi ^1' 1' 1' ^4; 2:41, 1, 1, 1) X 

Zf^iix, y; 1, Z2, Z3, 1; 1, 1, 223, 1) 

— ^41 ^4]^ ^23 2:23 Z^_^{x, y; Zl, Z2, 1, 1; 1, 2^12, 1, 1) X 

^^!°i(a;, y; 1, 1, Z3, z^; 1, 1, 1, ^34)), (133) 

where x, y and zi, . . . ,Z4 are arbitrary, while the remaining parameters are given by (130)- 
(131). 

Finally, the required identity (79) follows from (133) by applying (104), including the 
special cases (105), to each of the six cases of generalized ASM generating functions, and 
using (130)-(131) to eliminate ^41, 2:12, -223, -234, ^41, -212, -223 and 2:34. 

5.7. Alternative statement of some results of Section 4.1. In this section, it is shown 
that some of the results of Section 4.1 can be stated, slightly more compactly, in terms 
of a certain function which will be defined in (134) (and has essentially already appeared 
in 133). Also, an additional symmetry property for ASM generating functions is obtained, 
and further quadratic relations satisfied by these functions are discussed briefly. 
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Define a function 

AeASM(n) 

Pt(A)+Aii-1 pr(A)+A„„-1 pb(A)+A„„-1 pl(A)+Aii-1 
^1 ^2 ^3 ^4 ^ 

(|/Z4^1-(^4-l)(^l-l))'"^" (y-x(^l-l)(^2-l))'"^"' X 

(2/^2^3-(^2-l)(^3-l))'"^"" (y-x(^3-l)(^4-l))'"^"^ (134) 

It follows that Yn{x,y; zi, Z2, Z3, Z4) is a polynomial in x, y and 2:1,..., 2:4, with integer 
coefficients. 
Also define 

Yn{x, y; zi, Z2) = F„(x, y; zi, Z2, 1, 1). (135) 

It can be seen, using the definition (97) of the generalized ASM generating function, that 

F„(x, y; zi, Z2, z^, Z4) = y^Zf °(x, y; zi, Z2, Z3, z^, 

]_ {Z4,-I){zx-1) -^ x{zi-l)(z2-l) -^ {Z2-1){Z3-1) ^ X{Z3-1){ZA-1) \ (']_3g) 



i.e., Fn(x, y; zi, 2:2, 2:3, z^) = j/^Zf "(x, y; zi, Z2, Z3, Z4; Z41, zu, Z23, ^34), with 2:41, 2:12, 2:23 and 2:34 
given by (130). 

Therefore, by applying (104) to the RHS of (136), Yn{x, y; zi, Z2, z^, Z4) can be expressed in 
terms of quadruply-refined, adjacent-boundary doubly-refined and unrefined ASM generating 
functions. 

It can be seen immediately from (134)-(135) that 

Yn{x, y; zi, 1, 2:2, 1) = y^ ^7^(x, y; zi, 2:2), 
F„(x,|/;z, 1) = y^Z„(x,?/;z), 
y„(x,y;l,l) = |/^Z„(x,y), (137) 

and by acting on ASM(n) in (134) with transposition or anticlockwise quarter-turn rotation 
and using (6) (or by applying (100) to (136)), it follows that 

F„(x, y] zi, 2:2, Z3, Z4) = Yn{x, y; 2:4, 2:3, 2:2, zi) 

= x"("-)/^+^ iz,Z2Z,z,r^ r„(i, f ; ^, ^, ^, ^). (138) 

It can also be shown, using (134) and the properties of ASMs in which a 1 on a boundary 
is in a corner or separated from a corner by a single zero (as discussed in Section 2.1), that 

F„(a;,7/;zi,Z2,^3,0) = (l + ^^)(l + ^^)rn-i(x,2/;zi,Z2,^3,l). (139) 
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Taking x, j/, zi,...,Zi, 2:41, 2:12, -^23 and 2:34 to be parameterized by (126), and us- 
ing (120), (130) and (136), gives 

b{r, s)("-i)("-4)+8 c{r, s)"-« (6(ti, s) b{r, h) h{t^, s) b{r, U)^ t}/' 4^' tj^' t^' (rs)-i x 

-i^ I / a(r,s)\'^ /'c(r,s)\2_ a(ti,s) 6(r,s) a(r,i2) fe{r,s) a(t3,s) fe(r,s) a{r,t4) b(r,s) 
" \^\ fe(r,s) / ' V b(r,s) / ' a(r,s) 6(ti,s) ' a(r,s) 6(r,i2) ' a(r,s) b(t3,s) ' a(r,s) fe{r,t4) ^ 

Z^V(t,,r,...,r,t3;t4,s,...,s,t2), (140) 

where the weights in the partition function are given by (114), and the functions a, b and c 
are given by (121). 

It now follows from the symmetry of the partition function on the RHS of (140) in ti 
and ts, and in ^2 and ^4, as discussed after (123), that 

y„(x, y; zi, Z2, Z3, Z4) = Yn{x, y; Z3, Z2, zi, Z4) = F„(x, y; zi, Z4, Z3, 22). (141) 

In contrast to (138), there does not seem to be a simple combinatorial derivation of (141). 
(Note that the first equality of (138) together with either equality of (141) gives the other 
equality of (141).) 

It can be seen that (141), together with (104), provides an identity (referred to after (90)) 
involving Z^'^^'^{x, y; Zi, Z2, Z3, za) and Z^"^'^(x, y; 2:3, 2:2, zi, za) (or Z'^'^'^{x, y; Zi, Z2, Z3, za) and 

Z^"'^'i(a;,i/;zi, 2:4, 2:3,^2)). 

It follows from (133) and (136), together with (130)-(131), (135) and (137)-(138), that 

^(2:4-22) (2;i-2;3)y„(x,|/;2;i,22,2;3,2;4)Z„_2(x,|/) = 

(x(2;i - 1) (2:2-1)-!/) ((2:1-1) (2:2 -l)-i/2;i 2:2) (x(2;3-l) (2:4 -1)-|/) x 

((z3-l)(2:4-l)-y2;32;4)r„_i(a;,y;zi,Z4)F„^i(x,2/;2;2,2:3) - 
(a;(z4-l)(2;i-l)-y)((z4-l)(2;i-l)-y2;42;i)(x(z2-l)(2;3-l)-y) x 

{{z2-l){z2,-l)-yz2Z3) Yn-i{x, y, zi, Z2) Yn-i{x, y; Z3, za). (142) 

This can be regarded as a restatement of (79), where (104) and (136) enable the conversion 
between (79) and (142). 

Setting 2:2 = 1 in (142), and relabeling 23 as 2:2 and za as 2:3, gives 

y{.zi-z2) (2:3-1) ^n(a;, y, z^, 1, 2:2, 2:3) Z„„2(a;, y) = 
{x{z2-l)iz3-l)-y){{z2-l){z3-l)-yz2Z3)ziYn-iix,y; Zi, Z3) Zr,-i{x,y; Z2) - 

{x{z3-l){zi-l)-y){{z3-l){zi-l)-yz3Zi)z2Yn^i{x,y; Z2, Z3) Zn-iix,y; zi), (143) 

which can be regarded as a restatement of (81). 
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Setting 2:3 = in (142), using (139), replacing n by n + 1, and relabeling z^ as Z2 and z^ 
as 23 gives 

y'^{zi-Z2) Z'i F„(x, y\ zi, 1, Z2, ^3) ^n-i(a;, y) = 

(2:1 -l)(x(2;2-l) (2:3 -l)-y) ((2:2-1) (2:3 -l)-l/2;223) F„(x,|/;2;i,2;3)Z„_i(x,|/;2;2) - 
(2:2 -l)(x(2;i-l) (2:3-1)-!/) ((2:1-1) (2:3 -l)-|/2;i2;3) Yn{x,y\Z2,z-i) Zn-iix,y; zi), (144) 

which can be regarded as a restatement of (84). 

It can be seen that various other identities from Section 4 involving quadruply-refined, 
triply-refined or adjacent-boundary doubly-refined ASM generating functions could similarly 
be restated in terms of the functions (134)-(135). 

Finally, note that some further quadratic identities satisfied by ASM generating functions 
can be obtained from (142), or by combining certain variations of the Desnanot-Jacobi 
identity with the Izergin-Korepin formula, for certain assignments of the spectral parameters, 
and then using (140). 

An example of such an identity is 

(2^1-^2) (23-^4) Yn{x, y; zi, wi, Z2, W2) Yn{x, y; Z3, wi, Z4, W2) - 

(zi-zs) (2:2-2:4) Yn{x, y; zi, wi, 2:3, W2) F„(x, y; 2:2, wi, 2:4, W2) + 

(2:1-2:4) {z2-Z5)Yn{x,y;Zi,Wi,Zi,W2)Yn{x,y]Z2,Wi,Z3,W2) = 0. (145) 

This reduces to (60) for wi = W2 = 1, and can be obtained by applying (142) to each of 
the six cases of {zi — Zj)Yn{x, y; Zi, wi, Zj, W2) on the LHS, and then checking that the overall 
expression on the LHS vanishes. 

5.8. Derivations of (19) and (20). In this section, the explicit expressions (19) and (20) 
for the ASM generating functions at y = x + 1 are derived using a method different from 
that used in Section 3.2, which instead involves the six-vertex model with DWBC 

Let the crossing parameter q, for the weights (121), be given by g = ±i, this corresponding 
to the so-called free fermion case of the six- vertex model. It can be shown that, for this 
assignment of q, the partition function (113), with weights given by (114) and (121), is 
explicitly 



^n (""l) • • • '""nj Un+l, ■ ■ ■ ,U2n)\g^_^^ 



nr=l u/^^y^ Ul<i<j<ni'^i + Uj){Vi + Vj). (146) 



This result can be obtained by combining the Izergin-Korepin formula (123) (in which q = ±i 
gives a{u, v)b{u, v) = —{u^ + f ^)) with the Cauchy double alternant evaluation 

1 ^ ni<i<j<n("i-«i)(A-/3i) 



det ^—- = ^^^^'^^"^ , ^^V^ • (147) 
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For previous appearances of (146), see, for example, Bogoliubov, Pronko and Zvonarev [20, 
Eq. (58)], or Okada [Hi, Thm. 2.4(1), third eq.], and for information regarding (147), see, 
for example, Muir [ ■ , p. 345], [109, Sec. 353]. 

Assigning the spectral parameters in (146) according to (115), applying the relation (120), 
and using the first ten paramaterizations of (126), which enable arbitrary x and zi, . . . , 2:4 to 
be obtained, with y, Z41, Z12, Z23 and 2:34 then being given hj y = x + 1 and (130), it follows 
that 

7gen /'™ ™ I 1 . ^ ^ ^ ^ . 1 (Z4-I)(zi-1) i x(zi-l)iz2-l) i (z2-l)(z3-l) i x{z3-l){z4-l) \ 
Z„ [X,X+ L,Zi,Z2,Z3,Z4,L- (^+i) ^^ ^^ ,i J^ '^^ {x+1) Z2 Z3 ' "^ I+l ) 

= {xZiZ3 + l){xZ2Z4 + l)[{xZi + l){xZ2 + l){xZ3 + l){xZ4 + l)y X 

(a;+l)(n-2)(n-3)/2-2n+3_ (^43) 

(Using (136), the function in (148) could be written as {x + l)~^Yn{x, a; + 1; zi, Z2, z^, Z4).) 

The required expressions (19)-(20) now follow from (148) by applying (104)-(105), and, 
in some cases, setting boundary parameters to 1. 

5.9. Derivation of (71). In this section, a derivation is given of the identity (71) satisfied 
by the ASM generating function (66) associated with several rows (or several columns), for 
certain assignments of some of its parameters. 

For < m < n and 1 < ki < . . . < km < n, and indeterminates r, s and ti, . . . ,tm, 

consider the partition function Z^{r, . . . ,r,ti,r, , r, tm, r, . . . ,r; s, , s), as given 

by (113)-(114), with arbitrary functions a, b and c, where tj appears in position ki within 

r, . . . , r, ti, r, , r, tm, r, . . . ,r. Using the bijection of Section 5.2 between ASM(?2) and 

6VDW(n), and the behavior (111) of the relevant statistics under this bijection, it can be 
checked (by applying a process similar to that used in Section 5.2 for the derivation of (120)) 
that this partition function is related to the ASM generating function of (66) by 

^n y-i ■ ■ ■ i^i'l^ii^i 1^1 tm, f,...,r;s, , s) = 

b{r, s)(""i)(«-'") c{r, s)"-'" n^^ b{ti, s)"-i c{ti, s) x 

/7fci,...,fc^ f/ a{r,s) \2 / e(r,s) \2_ a{ti,s)b{r,s) a{tm,s)b{r,s) _ / c(fi,s) \2 / c{tm,s) \'2\ ( T AQ\ 

^n \\h{T,s)J ^\b(r,s)) ' a{r,s)6(ti,s)' • • • ' a{r,s)fe{t„,s)' U(ti,s)/ ' ' ' ' ' U(im,s) / J' \^^^ J 

(Note that the case of (120) obtained by setting ^2 = ^4 = s, and then relabeling ^3 as t2, 
matches the case m = 2, ki = 1 and k2 = n oi (149).) 

It follows from a result of Colomo and Pronko [47, Eq. (6.8)], [iS, Eq. (A. 13)] (with a 
special case stated previously by Colomo and Pronko [46, Eq. 5.8]) that if the functions a, b 
and c are given by (121) (which involves the further indeterminate q), then, for 1 < m <n, 



(a(r,s)6(r,s))™(™-i)/2 ni<.<,<,n(^^ - ^.) x 

Z^V(r, . . . , r, ti, r, , r, t^, r, . . . , r; s, , s) HIl'I^ Z^-ii^, . . . , r; s, . . . , s) 

= det (ia{t„s)biU,s)y~'{U-rr~^Zll^^,iU,r,...,r;s, ,s)), (150) 
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where the positions of ti, . . . , t^ within r, . . . ,r,ti,r, , r, tm, r, . . . , r on the LHS of (150) 

are now immaterial. An alternative proof of (150), in which a general identity for minors 
of a matrix is combined with the Izergin-Korepin formula (123), will be given elsewhere by 
this author. A result related to (150), in the case q = e^^'^*/^, has been obtained by Ayyer 
and Romik [12, Thm. 5]. 

Using (149) to express each case of a partition function in (150) as a case of the ASM 
generating function (66), and applying the definition (69), gives 

ryki,...,km (( a(r,s) \'i / c(r,s) \2_ gfti, s) fe(r,s) a{tm,s)h{r,s) _ / c(fi,s) \2 / c(tni,s) \2\ _ 

^n \yb{r,s)) ^\b(r,s)) ' a{r,s) 6(ii,s) ' ' " " ' a{r,s) b(i„,s) ' \b(ti,s)) ' ' ' ' ' U(i™^,s) J ) 

Y f(2i££l\'^ / c(r,s) \2_ a{ti,s)b{r,s) a{tm,s)b{r,s) \ /i ci \ 

^■^\\b(r,s)) ^\b{r,s)) ' a(r,s)6(ii,s)' • • • ' a{r,s)fe(i„,s)y • K^'^^) 

The required result (71) follows immediately from (151) by parameterizing x, y, zi, . . . ,Zm 
and Wi,...,Wm in terms of q, r, s and ti, . . . , t^ as x = (fM) , y = (^M) "• " °'^*"'''' ^^'"'''' 



b{r,s) J ' ^ \b{r,s)J ' * a{r,s)b{ti,s) 



and Wj = ( w^"^| ) , and observing that this parameterization enables arbitrary x, y and 
2;i, . . . , 2;^ to be obtained, with wi, . . . , Wm then being given by Wi = xzf + {y — x — l)zi + l. 



5.10. Derivations of (48) and (73). In this section, an expression for the partition function 
of the six- vertex model with DWBC at a certain value of its crossing parameter is given. This 
result, together with certain other previously-stated results, is then used to obtain derivations 
of the identity (48) satisfied by the alternative quadruply-refined ASM generating function 
at X = y = 1, and of the identity (73) satisfied by the function (69) aX x = y = 1. 

Let the crossing parameter g, for the weights (121), be given by g = e^^'^*/^, this corre- 
sponding to the so-called combinatorial point of the six- vertex model. It can be shown that, 
for this assignment of g, the partition function (113), with weights given by (114) and (121), 
can be expressed as 

(±1)"^"' S^/^ul^^ ..Uli' S(„_i,„_i,...,2,2,l,l)(Wl, . . . ,M2„), (152) 

where S(„_i,„_i_...,2,2,i,i)(^i, • • • ,M2n) is the Schur function indexed by the double-staircase 
partition (n — 1, ra — 1, . . . , 2, 2, 1, 1), evaluated at the spectral parameters ui, . . . , U2n- 

An important consequence of (152) is that Z^ {ui, . . . , Un] Wn+i, • • • , W2n)|g=e±2Ti/3 is sym- 
metric in all spectral parameters Mi, . . . , M2„. 

The result (152) was first obtained by Okada [ , Thm. 2.4(1), second equation], and 
Stroganov [: :, Eq. (17)]. For further related information, derivations and results, see, for 
example. Aval [ ], Fonseca and Zinn- Justin [75, Sec. B.2], Lascoux [ ;8, p. 4], Razumov and 
Stroganov [117, Sec. 2], or Zinn- Justin [149, Sec. 2.5.6]. 
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Proceeding to the derivation of (48), using the symmetry of the LHS of (152) in all of its 
spectral parameters, it follows that 



Z„ (ti, 1, . . . , 1,^3; *4, 1, • • • , 1,*2)| 



q=e 



±2Tvi/3 



Zn (^1,^2,^3;^4, 1, • • • , 1; 1, ! 1) L_p±2,ri/3 • (1^3) 



Applying (120) and (149) (with r = s = l, m = 4 and ki = i) to the LHS and RHS, 
respectively, of (153), and then using (121) (with q = e='=^'^*/^) for the functions a, b and c 
(and noting that a(l, 1) = 6(1, 1) and 6(1, 1)^ = c(l, 1)^), gives 



/ fe(l,t2)b(l,t4) \"-l/7gen/-| i. a(ti,l 
\b{t2,l)b{U,l)) ^n \^-L5 -L, h{ti,l 



a(ti,l) a(l,t2) a(t3,l) ajlM) . 
)' b(l,i2)' 6(i3,l)' b{l,t4)' 

a{l,l)a(ti,ti) b{l,l)b(ti,t2) a(l,l) a(f3.t2) b(l ,1) bjta fy) 

a(ti,l)a(l,t4)' b{ti,l)b{l,t2)' a(i3,l)a{l,i2)' 6(t3,l) fe{l,t4) 



Q^gzh27ri/3 



71,2,3,4/1 -,. a(ti,l) a(t2,l) a(f3,l) a(t4,l) . / c(tl,l) \2 /^ c(f2,l) \2 /^ c(t3,l) ^2 / c{t4,l) \A /. p- .\ 

^n \^-L) J-, b(ti,l)? 6(t2,l)' b(t3,l)' 6(t4,l)' Vb(ti,l)J ' U(t2,l)^ ' U(i3,l)/ ' V6(t4,l)^ J ±2^3' ^ ^ 



It follows from (154), by applying (151) to the RHS, and then setting U = ^^, or 
equivalently Zi = ^^^ = 7:,'",{ , for each 1 < i < 4, that 



(j^gzb27ri/3 
I 



(z2-24)""'^r( 1, 1; 21, -, Z3, ^; 1+ (fi^llkl^, ^ (£i^l){£2^ ^ (.2-l)(.3-l) l+ (^3-l)(.4-l) 

\ ^ ^/ " \ ■ ' -^' 22 ' ' 24 ' Zl ' 22 ' 23 ' 24 

= Xn{l,l;Zi,Z2,Z3,Z4^). (155) 

The required result (48) now follows from (155) by applying (104) to the LHS, and using 
the first definition of (10), the definition (69) and the identity (89). 

Proceeding to the derivation of (73), setting r = s = 1, q = e^^'^*/^ and k^ = iin (150), and 
then applying (152) to the partition function Z^{ti, . . . , t^n, 1, • • • , 1; 1, • • • , 1), applying (149) 
(with m = or m = 1) to the remaining cases of partition functions, and using (69) and (121), 
gives 

S(n-l,n-l,...,2,2,l,l)(^l, • • • , ^m, I5 • • • 5 1 ) = 

2n—m 

3n(n-i)/2 ^(1^ i)-„.(n-i) (^(^^^ ^y^^ ^(^^^ i))„-i ^^(^^ 1; f|;i}, . . . , fl^) |^^^,2W3- (156) 

(Alternatively, this can be obtained by combining (121), (149), (151) and (152).) 

The required result (73) now follows from (156) by setting tj = ^x^, or equivalently 

_ qti-l _ "U. ) -fQj, gg^(.]^ I < i < 771. 

* q-U b{ti,l) ' — — 
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